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! 

Studies were conducted on the large deformation response of 
composite beams subjected to a dynamic axial load. The beams were 
loaded with a moderate eccentricity to promote bending. The study was 
primarily experimental but some finite element results were obtained. 

* Both the deformation and the failure of the beams were of Interest * The 

static response of the beams was also studied to determine potential 
differences between the static and dynamic failure. Twelve different 
laminate types were tested. The beams tested were 23 in. by 2 in. and 
generally 30 plies thick. The beams were loaded dynamically with a 

gravity-driven impactor traveling at 19.6 ft. /sec. and quasi-static 
tests were conducted on identical beams in a displacement controlled 
manner. For laminates of practical interest, the failure modes under 
static and dynamic loadings were identical. Failure in most of the 
laminate types occurred in a single event involving 4QX to 50£ of the 
plies. However, failure in laminates with 30 c or 15’ off-axis plies 
occurred in several events. All laminates exhibited bimoduiar elastic 
properties. The compressive flexural moduli in some laminates was 
measured to be 1/2 the tensile flexural modulus. Mo simple relationship 
could be found among the measured ultimate failure strains of the 
different laminate types. Using empirically determined flexural 
properties, a finite element analysis was reasonably accurate in 
predicting the static and dynamic deformation response. 
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Chapter 1 
INTRODUCTION 


Fiber reinforced plastic composite materials have seen increased 
use In the aircraft industry in the past decade. These materials have 
found many applications in secondary and primary structures of military 
aircraft. Currently their use on commercial transport aircraft is 
limited to secondary structures. However, the use of composites in the 
primary structural design can significantly reduce the weight and 
improve the fuel efficiency of an aircraft. So it is likely that 
composite materials will see increase usage in the primary structural 
design of future commercial transport aircraft. 

Much of the research in composite materials has been in determining 
the structural response ard fatigue life under in-flight conditions. 
With commercial air transport an additional concern must be addressed. 
This additional concern is crashworthiness. Crashworthiness is 
concerned with preserving the well-being of the crew and passengers 
during crashes. The crashworthiness of a structure involves many 
issues. Fuel containment, seat design, peak deceleration, preservation 
of occupant space, body motion restraint, flammability, and smoke 
toxicity are just a few (ref. 1). However the main requirements of a 
crashworthiness structure are to maintain a protective shell for the 
occupants and to reduce their peak decelerations (ref. 2). Since 
crashworthiness deals with potentially survivab'e crashes, the impact 
velocities considered are relatively sna'-i (less than 40 ft/s). With 
impact velocities greater than this, the chances of surviving are 


1 


2 


minimal . 

Because of elastic-plastic behavior, aluminum absorbs energy when 
it yields and fails. On the other hand, fiber reinforced plastics fail 
in a brittle manner, thus, in general, reinforced plastics have less 
capacity to absorb energy. However, the total energy absorbed by the 
structure is not necessarily the crucial factor for crashworthiness. 
Crashworthiness is concerned with designing failure processes that 
maintain a constant load as a function of deformation during t-he impact 
event. The constant load minimizes decelerations, yet absorbs energy in 
a steady and progressive fashion (ref. 1). Much of the crashworthiness 
research of composite materials to date has been or. testing the energy 
absorption capabilities of axially crushed tubes and specially designed 
honeycomb structures. It has been found that by changing the fiber 
orientations, material systems, and the structural gentry, the 
stability of the collapse can be controlled. Consequently the energy 
absorption can be significantly increased with a stable collapse. These 
tests give useful Indications of the energy absorb ‘ng capacity for 
crushable parts of a structure. However, it is not entirely clear at 
this point in the development of composite structures that special 
crushable structures are necessary or all aircraft to meet 
crashworthiness criteria. Before designing special crushable structures 
It is necessary to know now a composite fuselage, designed for 
structural efficiency under in-flight conditions, will respond to crash 
loadings. Will the deceleration of the occupants be severe enough to 
cause death or injury? Will the integrity of the seat area be 
maintained? will fuel leakage be a isjor problem? Khai, in general. 
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win be the dynamic response of a composite aircraft to crash 
conditions? Of course the dynamic response of a fuselage will vary 
greatly, depending on the attitude and the velocity of the airplane upon 
impact, and by changing the structural configuration of the fuselage. 
One way of understanding and predicting the response of the fuselage 
under the many ways a survivable crash could occur would be to perform 
large deformation dynamic aralysis using a nonlinear finite element 
program. However, little is known about the load-deformation behavior 
of composite materials under dynamic loading and large deformations. 
Therefore, some of the basic input information for such an analysis is 
not available. 

Very little researcn has been done on the large deformation 
response and failure of laminated composite structural elements, whether 
static or dynamic, let alone full composite structural assemblies. 
Understanding the response of simple structural elements is essential to 
being able to perform successful crash analyses of a complete composite 
fuselage. It is certainly the least expensive way to approach the 
problem, both experimentally and computationally. 

Full scale aluminum fuselage sections of transport aircraft (Boeing 
707) have been vertically drop-tested (refs. 3 and 4), The correlation 
of test data and the finite-element model was quite good for the global 
deformations and the decelerations. From the results of these studies, 
it appears that the large deformation dynamic bending response of 
aluminum structural elements Is well understood. In addition, the 
vertical drop tests show that the majority of energy absorbed during a 
crash test is dee to bending failure of skin, stiffeners and stringers 
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(ref. 5). Little or no energy is absorbed by the stable progressive 
crushing of components. 

In this spirit, then, this study investigated the large deformation 
bending response of composite beams with a rectangular cross section. 
These beams could have represented the fuselage sScin, cr stiffener 
flange, or stiffener cap. The beans were loaded dynamically in bending 
to simulate the situations observed in the crash testing of aluminum 
fuselage sections. This study is considered to be a first step in 
understanding the dynamic response of composite structures. The overal' 
goals of the study were to: 


1. Design a simple test fixture to introduce crash-related bending 
loads in beams. 

2 . Determine any difference between large deformation static and large 
deformation dynamic response. Specific interest is in the failure 
mode. 

3. Determine the influence of laminate stacking arrangements on the 
dynamic response and failure mode. 

4. Predict the static and dynamic response, using an existing finite 
element program. 


The majority of the effort in this study was experimental. To 
successfully simulate large deformations and failure under crash related 
loads, a suitable test fixture was designed and built. The fixture was 
built around the concept of a drop tower. The fixture is discussed in 
Ch. 2 . Next, to determine the effects of dynamics on the large 
deformation response, the static resper se had tc be understood, or at 
least observed. So, re any dynamic testing was done, three beams of 
a given laminate type were tested under quasi-static loads. The load, 
deflection, strain, and failure response were observed. 
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To determine the effect of laminate type on the response* twelve 
different laminate types were tested. After the static tests were 
completed, three more beams of each laminate type were tested 
dynamically. The results of the static and dynamic experiments are 
presented in Ch. 3. The differences in the dynamic and static tests, 
and the differences In the response and failure modes between the 

different laminate types are reported. An unexpected result in the 
surface strain response of the laminates was observed. Accordingly 

additional tests were conducted. In Ch. 4, this test procedure is 

presented and the results are reviewed to verify and further understand 
the phenomena. 

Methods for predicting the response of the beams are explored in 
Chs. 5 and 6. Predicting the response of the beam involves both 

predicting the on-set of failure and predicting the global deformations 
to the applied load history. The strains at first failure in each 
laminate and the success of a strain-related failure criteria at 
predicting these failures are examined in Ch. 5. An existing finite 

element program was used in Ch. 6 to determine the success of predicting 
the deformation response of the composite beasrs. The predicted time 
histories of force and displacement response are presented and compared 
with experiment, as are the static and dynamic load-displacement 
response. In addition the spatial shapes of the dynamic response is 
presented at various times after impact. Finally, Ch. 7 presents 
conclusions and recommendations for further study. 


Chapter 2 

EXPERIMENTAL SET-UP ANO PROCEDURE 


The basic leading configuration used in the study is shown ir fig* 
2 . 1 . The beam was loaded In a column fashion with a moderate amount of 
eccentricity to promote bending. The beams were oriented vertically as 
shown. The lower end of the beam did not move while the upper end moved 


vertically. A hinge on each end allowed free rotation of the beam 
ends. The static tests were conducted in a displacement controlled 
manner in a standard screw driven load frees. The dynamic tssts were 
conducted in a specially designed drop tower and the dynamic loading, 
denoted by F(t) in fig. 2.1, was provided by a gravity criven 
impactor. In the unloaded position the hinge supports were 24.2 Ir. 
apart. The specimens were 23 in. long by 2 in. wide and generally 30 
plies thick. At each end 1.5 in. of the specimen was clamped In the 
hinge, leaving the unsupported length of the t8S" 2t 20 in. In tctfi the 
static and dynamic tests the deflection of the upper end of the beam was 
limited to 16 in. by a bumper. This limit on deflection was iirpcsec so 
there would not be damage to the test fixture or related instrumentation 
due to the metal hinges impacting each other. 

A 4 point bending configuration was originally considered. 
However, the eccentrically loaded column configuration *as used 
instead. The primary reason was that to successfully model crash 
conditions, large deformations must be irtroduced into the beams. For 
large lateral deformations there are considerable axial deformations as 
well. With a 4-poirt loading system, the beam could easily deform to 
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such an extent that it would slip through the supports without failure, 
as shown in fig. 2.2. With the eccentrically loaded column the beam 
deflections were limited only by the ends of the beam touching. In 
addition, with the 4-point loading systems ths bssR srsds sisovs rs! stive 
to the simple supports. As the beam moves over the support, 
considerable f riction forces can be generated. Tnese forces are 
difficult tc model analytically. Cn the ether hand the eccentrically 
loaded column's simple support moves with the beam end. In this case 
the support can be attached to a low-friction linear bearing. 

Another draw back of the 4-point loading system is that it can 
cause failure at the supports. With a four point bending configuration 
the supports impact on the surface of the beam in the regions of maximum 
bending moment. The impact can cause local damage and Initiate 
failure. With the eccentrically loaded column the two support points 
are on the ends of the beams, regions of minimum bending moment. 


Consequently, failure always initiates at the center of the beam away 
from the end. With failure occurring in the center, the complicating 
effects of stress concentrations at the supports are not present. 


2.1 Beam Specimens 

Table 2.1 shows the lay-up and number of plies for each of the 
laminate types tested. The laminates were fabricated by the JJASA- 
Langley Research Center usirg AS4/3502 graphite -epoxy pre-prea tape. 
For each laminate a 24 x 24 in. panel was fabricated. After curing, the 
panels were C-scanned to determine if any defects were present. Then, 
ten 2 x 23 in. seams were cut from each panel. Throughout tnis report 
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Beam Slipping 


Figure 2.2 A Possible 4-Point Loading Arrangement 
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TABLE 2.! 

SPECIMEN LAMINATE ORIENTATION 


LAMINATE 

LAY-UP 

NUMBER OF 

IMPACTOR 

NUMBER 


PLIES 

WEIGHT 

(LBS) 


I 

t0| 30 

30 

48.6 

2 

l ( 15/0/-15) 5 ) s 

30 

36.1 

3 

I ( 30/0/— 30) 5 1 s 

30 

24.5 

4 

[(45/0/-45) 5 J s 

30 

24.5 

5 

[ ( 60/0/-60) 5 1 s 

30 

24.5 

6 

[ ( 75/0/-75) 5 1 s 

30 

24.5 

7 

{(90/0/-90) 5 ) s 

30 

24.5 

8 

{ (0/90) 8 1 s 

32 

36.1 

9 


32 

36.1 

10 

U90 8 /0 8 ] $ 

32 

36.1 

11 

[ (45/-45/0/90) 4 l s 

32 

24.5 

12 

I(0/45/0/-45) 3 /90/0/0 1/2 I s 

29 

36.1 



n 


etch beam specimen will be referred tc by its designated number, e.g., 
specimen 2.1. The whole number portion refers to the laminate type, and 
the decimal portion refers to the particular beam cut from the panel. 
For example, specimen 2.1 ts beam no. 1 cut frcsss the { (15/0/-15) 
panel. 


The laminate types were separated into three groups for study. In 


the first group, numbers 1 tnroug 


mttnk 7 


k Ka *1 

• iiC iQiHmubCa 


.iAwa * 1 

wet c u » ; wi 


the ((e/0/-9) 5 i s family. The angle 9 increased from 0° to 90“ In 

steps of 15‘ to study the effect of increasing the angle of the off- 
axis laminae. The second group, "uSbei’S 8, 9 and 10, wsrs orthctrcpic 
laminates which all had the same inplane stiffness. The plies were 
stacked differently to study the effect of clustering on the failure 
modes and response. In the final group, nunbers 11 snd iwSrs 
laminates which, according to current design philosophies, might be 
found on a fuselage. Laminate 11 is a quasi -isotropic lay-up and 
laminate 12 had 52 % 0's, 41% 45' s and 7* 90 s. Lass' nate 12 is issue h sers 
orthrotropic than laminate 11. 


2.2 Dynamic Test Fixture 

Figure 2.3 shows a drawing of the drop tower fixture used for the 
dynamic testing. Two 10 ft. long vertically oriented hardened steel 
rods 1 in. in diameter spread 6 'n. apart were fastened tc C channels ers 
either end. The channels were f’xec to the floor and celling of the 
room used for testing. The slide*", which suoported the uppe*" end of the 
composite beam, and the mass car, which provided the impact force. Slid 


down the rods on low-friction linear bearings. The composite beam 
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specimen was clamped to the tcp and bottom hinges which provided the 


a.y *, r\r% f U ^ 

un i ^ V* U'lC 


simple supports. Figure 2.4 shc*s a detail dr 
hinge. As shown in fig. 2.4, the rear surface of the specimen was 
offset 5/S *n. fr on- the hinge pivot, or pin, which was in the plane of 
tne verf'cai hardened steel reds. Th.*s offset provided the eccshtr ! city 
which promoted the bending, "ne tcc hinge was attached to the slide** 
whereas the bottom hinge, which was identical to the top one, was 
attached to the lower c charnel. Weights were fastened to the ?.ass 
car. The car was raised 6 ft. above the spherical steel impact point on 
the slider. A solenoid-activated release mechanism released the mass 
car. The mass car impacted the slider at 19.6 ft/s. As the sl ; d£- and 
the mass car combination move downward, the specimen deflected axially 
and laterally. Figures 2.5 arc 2.6 are pnotographs of the apparatus and 
instrumentation. 


2. 3 Static Test Fixture 

To conduct the static tests, the tcp and betters hinges from the 
dynamic test apparatus were attached to the heads of an lrstrcn load 
frame. Figure 2.7 shows specimen 2.1 partially deflected and failed in 
tne static test fixture. Aga’n, the beam v?as clanged to the hinges and 
offset from the pivot points by 5/8 *n. As can oe seen in fig. 2.7, tne 
hinges rotated with the enas or tne oeam. In tne 'oad frame the too 
hinge was attached to the "'cad ce 1 ! and r.as stationary. The botetm 
hinge was attached to the moving crcsshead. 
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Force 

Transducer 


5 /" 



Force 

Transducer 



Beam Specimen 


Figure 2.4 Sc*e»atic of Upper Hinge 









gure 2.7 Photograph of Static Test Fixture 
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2.4 Data Aquisition 

For both the static and dynamic tests, end load, end displacement, 
and surface strains at the center of the bean were recorded. The 
strains measured bending and axial compressive effects in the beam. 
Poisson strains were not measured. The manner in which these 
measurements were recorded differed in the static and dynastic tests. In 
addition to these measurements, for the dynamic tests high speed motion 
pictures of the edge view of the laminate were taken to record the 
overall deformations and the failure scries . 

2.2.1 Static Tests 

in the static tests the end axial load {F(t> in fig. 2.1) was 
measured using a 1000 pound load cell. The end axial displacement was 
measured and controlled with the crosshead speed which was set at 5 in. 
per min. The load - end displacement re laticn was recorded using the 
chart recorder. Back-to-back strain gauges measured the surface strains 
at the center of the beam and the signals were conditioned using 
amplifiers designed and built at the Langley Research Center. Strain- 
displacement relations from the two gauges we-e recorded using a X- v 
plotter. All data were then digitized manually to facilitate data 
reduction. 


2.2.2 Dynamic Tests 

The dynamic tests were conducted at the NASA Langley Research 
Center's Impact Dynamics Research Facility. The oata aquisition system 
there was designed to permit the simultaneous recording of 90 data 
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channels on one 28 track magnetic tape recorder using a constant 
bandwith FM multiplexing technique. This experiment required the use of 
5 data channels. The channels were: twc loads, two strains, and one 
displacement. 

For the dynamic loading, the axial end load was measured by the two 
piezoelectric force transducers shown in fig. 2.4. The transducers were 
connected in parallel and placed on both ends of the hinge pin. This 
arrangement measured the sum of the slider forces acting on both ends of 
the hinge pin and therefore the total axia^ force time history at the 
upper end of the beam. An Identical set of transducers were placed on 
the bottom hinge. The signals from the transducers were conditioned by 
charge amplifiers. These amplifiers produced analog signals that were 
recorded by the FM system. Again, two strain gauges were placed back- 
to-back on the center of the beam specimen. These signals were 
conditioned by amplifiers designed and built by Langley Research Carter, 
but different from the ones used in the static tests. 

"he beam end displacement was measured by an optical displacement 
transducer. The optical device did net require ary mechanical linkage 
fastened to the slider. The transducer was designed and built by 
Langley Research Center. The transducer worked as shown in fig. 2.8. A 
16 in. long tapered bar was fastened to the slider. Tn? outside edge of 
the bar tapered from a 1 In. width at the top to a 2 In. width at the 
bottom. A laser beam, which was diffracted into a horizontal line of 
light, shined on the front of the bar. Directly behind the tapered bar, 
a one-inch long array of 1024 light sensitive diodes was mounted 
horizontally. When the specimen was in the fully uprlgnt undeflected 


Tapered 

Bar 



Diode 

Array 




Transducer 
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position, the bottom of the tapered bar fully shadowed the diode 
array. As the top of the specimen deflected downward, more and more of 
the diodes became exposed to the light. When the specimen traveled Us 
full 16 in. of axial displacement, the entire array was exposed to the 
laser light. An electronic circuit counted the number of diodes exposed 
to the laser light. The counter generated an analog signal proportional 
to the number cf the activated diodes. With proper calibration* the 
optical transducer provided a signal proportional to the vertical end 
displacement of the beam. 

As mentioned previously there were five t^me histories recorded 
during the dynamic event. These five dynamic signals were filtered at 
1-KHz to remove spurious noise. The filtered signals were connected to 
voltage-controlled oscillators where they were converted into descrete 
FM signals. The FM signals were then recorded on magnetic tape. After 
the testing was completed the magnetic tape was played back through the 
voltage-controlled oscillators and the resulting analog signal was 
digitized at 4000 samples per second. 

In addition to the digital data, a high speed movie cansera was used 
to record an edge view of the laminate as it deformed and failed during 
the dynamic event. The edge of the laminate was painted white to 
facilitate viewing. The 16mm camera was set up to take 4CC 
frames/sec. Movies were made for at least two soecimens from each 
laminate type. 

2.5 Initial Measurements 


Some initial measurements were taken of the geometry cf the beam 
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specimens. The results are shown In Table 2.2. From each laminate type 
the dimensions of 3 beam specimens were measured and then averaged. The 
length and width were measured to assure each was as specified. The 
thickness of each specimen was measured at the center and at eacn end of 
the beam with a micrometer. As can be seen, there was a variation In 
the thickness with lencth. The effect of the uncertainty In the 
thickness measurement on the response of the beasss is examined 1r 
Appendix A. The beams were next placed on a flat surface and the 
camber, or deviation from prefect straightness, was measured with a 
ruler to determine the initial eccentricities in the bean’s. The 
eccentricity would add or subtract to the built-in 5/8" eccentricity of 
the test fixture. However, the initial eccentricities in the beams 
could not be measured with absolute certainty. The effect cf the 
uncertainty in the eccentricity measurement is examined in Appendix B. 
Finally the beams were weighed. By knowing the Areal weight of the 
prepreg and the fina- weight cf the bean*, a good estimate of the fiber 
volume fraction could be made. As can be seen, the fiber voljve 
fractions were within expectations. It should be mentioned that no 
tests were conducted cr laminate 6. 



TABLE 2*? 

MEASURED SPECIMEN GEOMETRY 


LAW I HATE WIDTH 

HUMBER 


LENGTH 




AVERAGE 

THICKNESS 

ao.oo2* 
1TB 1ST 


CAMBER 

1 O.Ol* 
(Ip) 


FIBER 

VOLUME 

FRACTION 

X 


1 

2- 00 

23.00 

0.162 

2 

2.00 

23.00 

9.162 

3 

2.00 

23.00 

0.1SB 

4 

2.00 

23.00 

0.160 

5 

2.00 

23.00 

0.160 

6* 

2.00 

23.00 

— 

7 

2.00 

23.00 

0.164 

8 

2.00 

23.00 

0.174 

9 

2.00 

23.00 

0.175 

ID 

2.00 

23.00 

0.174 

11 

2. CO 

23.00 

0.178 

12 

2.00 

23.00 

0.154 


0.168 

0.162 

0.03 

69.6 

0.165 

0.162 

0.G2 

69.2 

0.160 

0.159 

C.02 

69.6 

0.165 

C.106 

0.02 

69.7 

0.167 

0.160 

0.04 

69.6 




— 

0.165 

0.163 

0.04 

69.8 

C.175 

0.175 

0.01 

68.2 

C.177 

0.176 

0.01 

68.2 

0.175 

0.174 

0.01 

68.5 

0.180 

0.178 

0.07 

69.1 

0.156 

0.155 

0.06 

71.2 


♦not tested 
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failure. Three were chosen siraply to determine the scatter *-= **** 


2.6 Test Procedures 
2.6.1 Static tests 

T hree specimens of each laminate panel were tested statically to 

v- 1 1« ♦!■*<"% <4 * ¥ i 
> tn wi c uava, 

particularly the failure data. The static test procedure was somewhat 
unique and is discussed here. first the load frame loadcell was 
calibrated. Next, the specinen was placed in the fixture such that tie 

initial camber in the beams increased tne load eccentricity. Then the 
load cell was zeroed, the strain gauges were connected to the 
amplifiers, and the amplifiers were balanced. The crosshead speed was 
set at 5 in./min. To determine if ro"visib!e damage like microcracking 
and fiber pullout was occurring as the beam specimens deformed, the 
axial end displacement was applied in stages, first the end of the beam 
was displaced axially 2 in. At the 2 in. displacement level, the 
cross head mot' on was reversed and the displacement of the beam end was 
returned to zero. Next the end of the beam was displaced axially 4 
in. At the 4 in. level tne crosshead motion was reversed and the end of 
the beam was returned to zero displacement. This procedure was repeated 


in 2 in. Increments until finally tne beam was displaced axially 16 in. 
and then returned to zero. If the loading and unloading load- 
displacement curves coincided, then, witnin the sensitivity of the 
Instruments, the specimen abso^bec no energy during the cycle. 
Conversely, if the two curves did not coincide, then tne area between 
the loading load-displacement curve and the unloading load- displacement 
curve was the energy absorfcec by the specimen due to the failure 
mecnanisms in that cycle. 
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2.6.2 Dynamic Tests 

Three more specimens of each panel were tested to failure 
dynamically after the static tests were complete. The 6 ft. drop he i gut 
selected for all the dynamic tests provided the same velocity of the 
mass car prior to impact for each test. However, since the stiffnesses 
of the 12 different laminate types were net the saiie, the decelerations 


of the mass at the top of the beam, and thus the dynamic force 'evels, 
would vary from one laminate type to the next. It was felt that it 
would be important for comparing the response of the different lasrinate 
types to have the dynamic force levels and the ve'ocities somewnat 
similar. To keep the velocity and force levels similar from one 


laminate type to the next, the ^inpactCT moss «ss changed or eacn 
laminate type. Tc detennine the impactor mass for each laminate type, 
the work required to deflect a beam 16 in. was calculatec from the 


static load-displacenert relation. An additional 2 G% rf energy was 
added which more than accounted for energy lost due to friction and 
during impact of the slide** and mass car. from this desired impact 
energy value, the impactor mass .vas computed. The impactc- rsess used 
for each laminate { s given in Table 1. 

Before any tests were attempted, calibration signals, corresponding 


to known levels of physical quantities, were recorded cn each char 


irh rhAnnpl trt 

^ , yv 


facilitate digitization. Viith the impactor mass in place, the specimen 
was clamped into the ninges and the 'instrumentation was connected to the 
recorders. Next, a 1 ! the transducers «ere balanced arc! zeroed. A 10 
sec. countdown procedure was used wherein at 5 sec. the FM cape 
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recorder was activated and at 1 sec. the movie earners was started. 
Finally, at time zero, the mass was released and the dynamic load, 
strain, and displacer it histories were recorded. 



Chapter 3 
RESULTS 


As described in the previous chapter, beam specimens of each 


laminate type were tested quasi-staticci ly an.j dynamical ly. The raw 
data from these tests were recjcec and put in a common format sc that 
easy comparisons of the test cond : tions arc the beam specimen •■esponse 
could be made. Comparisons were ticf between the static and dynamic 
tests *or each laminate type. Comparisons were also made aircrc lanrrate 
types. This chapter presents a discussion of these experimental resu'ts 
and related observations. 


3.1 Data Reduction 

3.1.1 Static Data 

'he raw static data were recorded J n graphic form, as mertionec 
previously, on an x-Y recorder. 'he graphic data were converted intc 
numerica' data by facing the cezves on 3 d ; citizirg tdC'c. T “Ch t"C 
digital oata was stewed ir computer f iles to facilitate potting and 
manipulation. tcad-displacemert and st**ain-cisp 'acement re'at'cns were 
generated from the cata to analyze t*e tests. 

3.1.2 Dynamic Jata 

The raw dynamic data were ’■ecc-roed on r M magret’C tace. T"? FM 
sigra’s were played back through vo' tage-control lec oscillators tc 
retrieve tne original analog signal. Ca‘*bratiom signals we^e used tc 
provide tne proper relations Detween the analog signals 2"d tne chjsica' 
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Quantities measured by the transducers. The analog signals were then 
converted to digital data at a sampling rate of 4QCQ points per 
second. The five digitized dynamic data channels from each test were 
then placed in a computer file. Upon initial scrutiny of all the data, 
both systematic and random errors in the displacement channel 5 were 
found. These data were enhanced to extract usable information as 
explained below. 

The random errors were believed to be associated with noise from 


the electronic circuit wh'ch produced the analog displacement signal. 
Figure 3.1 shows the displacement time history from specimen 1.4. Note 
that in the raw, or unenhancsd, signal, 8 < 2 m before the unpact occurred 
there was considerable noise. The no^se level In fig. 3.1 corresponds 
to axial disp'acement fluxuations of up to 0.4 in. With the impactor 
not in contact with the slider, these displacements are net possible. 
Therefore this portion o f the signal Is noise. Unfortunately, this 
noise persisted after tne impact and distortec the displacement 
signal. To improve this signal, the displacement data was passed 
through a digita' filter. The filter used was a nonrecursive lew-pass 
filte- with a cut-off frequency of 5C0 nz. A description of the filter 
des*gn and a **eview of its performance is given in Appendix u. 

In aedition to the random error, a considerable systematic error 
was ooserved as well. In fig. 3.1, in the raw signal, note the 

^frills^^nn 4 - A* irpirfimf hictnrw c AAn afro** t hp 

narmomc- 1 1 *ce csciiidLiOT ^ u ww * - * ,w 
impact point. An interpretation of this signal would indicate that in 
the first 0.01 sec of the event the beam end displaced downward 4 in. 
and then upward 0.5 1 r. This behavior would no: make physical sense. 
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The Inertia of the Impactor mass provided a monotonical 1y increasing 
displacement. Upward displacement was not possible. To check this, the 
high speed films were closely scrutinized. No evidence of the top of 
the beam traveling upward soor after Impact could be found in any of the 
films. The source of this error in the displacement signal is felt to 
lie in the displacement transducer design. Referring to fig. 2.8 and 
the tapered bar, in addition to up and down motion of the bar shadowing 
and exposing diodes, a small lateral displacement of the sl'der or the 


drop tower would also expose or shadow diodes. In fact, wth lateral 
motions of the slider the transducer would record an apparent axial 
displacement 16 times the lateral displacement. A small rotation of the 
slider in the plane of the drop tower would have a similar effect. 
There were several probable causes of this unwanted motion of the 
slider: First, if the stee’ spherical impact points of the slider and 
mass car were not exactly aligred, a lateral component of force or a 
moment would be imparted to the slider. This would cause lateral motion 
and prcbably rotation of the slider. This motion would be the result of 
tolerance in the bearings or actual flexural motions of the drop 
tower. Also, after the tapered bar of the displacement transducer was 
attached to the slider, the impact force no longer passed through the 
center of mass of the slider. This would definitely cause £ moment to 


be imparted to the slider. Unfortunately, during the design of the 
slider, the additional mass of the tapered bar was not taken into 
account . 

To remedy the systematic error, the first 2 in. of each 
displacement history was replaced by a straight me segment. Tne 
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segment was faired In by hand for each test. Additional ha**mon1c 
oscillations In the displacement curve after the failure point can be 
seen in fig. 3.1. These are cf such a frequency content that they are 
f e1t to be caused by the same unwanted motions of the tapered bar. 
However, these oscillations were not removed as they had no major 
effects on the results. The upper curve in fig. 3.1 shows the altered 
displacement time history used for specimen 1.4. The enhanced time 
history incorporated both the digital filtering and the Iritia 1 straight 
line segment. 

Once the displacement data has been enhanced, relations between the 
load, displacement, strain, and time data could be meaningfully 
plotted. Specifically, for each dynamic test the loac-displacement, 

strain-displacement, load-time, strain-time, and disp'acement-time 
relations are presented. 

3.2 Experimental Results 

To analyze each laminate type, seven relations «ere studied. These 
relations were: 

1) static load -displacement 

2) static strain-displacement 

3) dynamic load-displacement 

4) dynamic strain-displacement 

5) dynamic strain-time 

6) dynamic load-time 

7) dynamic displacement- time 
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In addition to these, the high speed movies were analyzed. The 
seven relations, the movies, and the post-test observations were used to 
characterize the test conditions and the beam response. The word 
'response' as used here means both the spatial deformations and the 
fal lure modes. 

By examining in detail the seven relations for laminate type 2, the 
( (15/0/-15) c l s laminates, the method of characterization of the test 
conditions and the beam responses will be described. Cue to the 
overwhelming amount of data from the tests, not all aspects of the 
testing of all the specimens will be described in detail. Laminate 2 
will be discussed ir detail and then a general description cf the 
results cc*won to all the tests wil 1 be given. Then a detailed 
description of the failure mode of each laminate type is presented. A 
comparison of the static and dynamic test results will be sade. 
Finally, a summary of the observations of all the testing will be 


provided. 

It should be noted, as previously stated, no tests were perfcr 
on laminate type 6. In addition, only two dynamic tests of laminate 
types 2, 3, 9, and 12 we-e completed because of equipment failure. 
Numerous strain guage failures were encountered on the static 
specimens. Typically only two static strain-displacement relations were 
recorded and in the case of laminates 5, 7, and 11 only one strain- 
displacement relation *»as recorded. 
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3.2.1 Experimental Characterization Method 

To characterize the static tests and the static beam response, 
relations between the end load and end displacement and relations 
between the two surface strains and end displacement were plotted for 
each of the static beam specimens. As explaired in the previous 
chapter, the displacements in the static tests were applied in 
increasing increments, and then returned to zero, to determine if 
nonvisible or nonaudible damage was occurring and dissipating energy. 
Figure 3.2 shows the load displacement relation for the specimen 2.2 
from the ( < 15/0/-15) 5 } s pane 1 . The horizontal axis shows the end 
displacement and the vertical axis the end load. During the first cycle 
the beam was displaced 2 in., then unloaded. This process traced and 
retraced the top-most curve. Since the unloading relation retraced the 
loading relation within the resolution of the instrumentation, it was 
assumed that no energy was absorbed in the beam during this first 2 in. 
of displacement. Again for the second loading and unloading, to 4 in. 
axial displacement, the relation traced and retraced the top-most 
curve. During the third cycle, to 6 in. displacement, a sha^p drop in 
the load occurred at about 5.8 in. of displacement. This sharp drop in 
the load corresponded to ply failures in che beam. These failures were 
definitely visible and audible. The ply failures reduced the flexural 
rigidity and thereby reduced the load required for that dlsolacement. 
The displacement was continued to 6 in. ano then unloaded back to zero 
load. With unloading, the load displacement relation followed the 
second curve from the top. As the loading-unloading cycles continued to 
increased lengths to 16 in., the loading-unloading curves coincided 


END LOAD. LBF 
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Figure 3.2 


Static Load-Displacement Relation 41th loading and 
Unloading Curves for a [ ( 15/0/-15) 5 l s laminate 
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except immediately after visible and audible laninate ply failures. 
This brittle failure-elastic response was typical for all the laminates 
except laminate 9. Therefore, it can be concluded that the dominate 
energy dissipation mechanism i n the beans was associated with -he Ply 
failure events. This is not to say that there was no nonvisible damage 
or that it did not dissipate energy. It is to say that energy 
dissipated by nonvisible damage could net 08 detected. do.*.s/er, ir tb- 
scale of energies associated with crash conditions, these nonvisible 
mechanisms absorbed negligible energy even if the mechanisms did 
exist. The unique static load-deflection behavior of laminate 9 is 
discussed later in section 3.2.3. 

For the sake of clarity and simplicity no other loading-unloading 
curves are presented in the main text. However, the loading'Unlcading 
curves for each static test are presented in Appendix D. Instead, only 
the outer most locus generated by the loading-unloading procedure is 
presented In the main text for each test. c iQjre 3.3 Shows OH 3 .omon 
plot the outer locus of the static ^ad-displacement relations for the 
three static specimens 2.1, 2 . 2 , and 2.3. The finite element prediction 
for the deformation response is also displayed cn fig. 3.3. th«jt 
prediction will be discussed in a later chapter. Note tne steep initial 
slope of the load-displacement relation as the stress state in the beam 
changed quickly from primarily colwsn-llke and compression to 
predominately bending. On the load-displacement relation, failure 
events are clearly indicated by sudden drops in the load. During a 
failure event several plies failed siinul tuMvOiiS^j cn 
the beam. The amplitude of the load drop gives a relative indication of 
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the number of plies which failed during the event. The similarity of 
the response of the three replicate specimens should be noted. 

Figure 3.3 also shows the static strain-displacerrent relations for 
the | ( 15/0/- 15) 5 1 s specimens. (Strain gage problems prevented recording 
two of the results from the 3rd specimen.) Since it was destroyed when 
the surface ply failed, the gage on the tension side of the bean 
recorded strain only up to the first failure event. The compression 
side guage continued to record. Progressive failures in the beam are 
clearly visible in these relations. The compress i «e strains suddenly 
drop because the ply failures reduce the thickness of the beam and 
thereby reduce the strain required for the curvature corresponding to 
that end displacement. ibe magnitude cf the strain levels shcud bs 
noted; strains in the 1-2* level were common. 

The high speed movies provided valuable visual Information on the 
dynamic response cf the beans. Figure 3.4 shows selected frames fret 
the film of soec’men 4.5, a ( (45/C/-45) 5 l s laminate. The approximate 
time, in seconds, after impact is indicated with each frame. The first 
frame shown, t = 0.0, was the last frame o* the fTm take r be f ore the 
Impact. Shortly after impact, the shape of the beam is quite 
interesting, as sncwn in the second frame of the film, approximately 
2.5 milliseconds after impact. In th*s frame the bears deemed intc 
a noticeable “W" shape. This shape occurs because the center of the 
beam has not yet responded to tne impulsive loading. The third frame of 
film shows the center of the beam as it s n a?s though th e 'u' 1 . The 
frame of film after that shows an acute curvature at the center of the 
beam. The final frame of film shews a more obtuse curvature in the 


Figure 3.4 Third Mode Deformation Response 
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center as the beam snaps back. This V shape and snapping Is 
associated with the third vibrat 1 or* ttiods cf ths bea*n» The ironies s> w« 
that In the test fixture the third mooe vibratory response in the beams 
was excited and that this response was superimposed on the global 
deformation response which was ■ ar tc tnc static rfsjjcrss. **«e ■ i 
of specimen 4.5 was chosen for f^g. 3.4 because it had the best visual 
information. However, the third node response was seen in all the 
dynamic tests. In fact, for a!: specimens the spatial defcrnjat'or.s up 
to first failure appeared similar. 

The upper portion of fig. 3.5 illustrates the load-displacement 
relations for the two dynamically testec [ { 15/0/- 15 ) ^ 1 s specimens, 
specifiers 2.4 and 2.5. The end displacement is shown on the horizontal 
axis and the end load on the vertical. The difference in scale Between 
the dynamic and static lead levels should be noted in. figs. 3.5 and fig 
3.3, respectively. The dyranre responses of the two beams were quite 
similar. The response curves up to failure a^e almost coincident, 
indicating tre high repeatability of the dynamic experiments. Due to 
the statistical nature o r ~ failure in corpcsites, there were some 
differences in the failure everts. As witn the static leading, failures 
are denoted by the sharp drops -n tho lead, Thprp an initial high 


amplitude load spike in the response as tne initial column configuration 
of the beam quickly decelerates the oeam. As the beam begins to bene, 
the spike subsides. The third mcce response, just discussed, is 
strongly evident in the lead response. However, this vibratory response 
dampened as the beam deformed axia’iy and laterally. This was evident 
in the f-iltrs, as well as in fig. 3.5. 


4 8 12 16 
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The lower portion of fig. 3.5 shows the surface strain-displacement 
relation for the two dynamic tests of the ( ( 15/0/-15>5 i s lamir * . *e 
convex (tension) side strain is on the tep ard the concave (compression) 
side stra ; n is on the bottom c c that portion of the f igure. *' „r mat 
the concave side (bottom curve) initially starts out in tension and tne 


convex side, (top curve) ' nitia j s.a^wj in ccmuress o •• 


"Trs’C '-niti,' 


reversal in strain is due tc the initial "*T deferred snape of tne 
oean. Tre fai'ure events A r e c'early evident on this figure. 

figure 3.6 shows tne load ard strain histories for the dynamic 
tests. The third node osculation can be further studied from this 
relation 'he frequency of oscillation can be determined directly from 
this figure and corparec to the computed 3^ node ritual frequency 
the beam. As w’th the previous figure, failure events a. e c'early 
evident. As ir the load- displacement re'ation, the load time history of 
the failure everts are denoted fcy sudde- crops in the 'ocC. K -C other 
•oad tine histories are presented *n me main text. The load tine 
histories of a'l the otner dyramic tests are presented in Append 1 ’ * f 

Tne displacement ti"e history for t"? two dynamic specimens tested 


is disp'ayed in fig. 3.’. Tre diso'aoenent ; s g : ven or vertical 
axis ard the t 'me cn tne norizerta’ ax's. Tre enhanced 2 °c raw cata a r e 
plotted fer 0 dc n s sre j ^ t s fre^ s*ch 520ri T £ n spdC0d 
vertically f rom eacn otne r in the figure. Note tne straight me 
relation in the enhanced oata rear time zerc. The end displacement-* 'me 
relation characterizes the cy n 2mic test. It provides inf ermat ’C r c r the 
total duration of the event and the approximate charge *n velocity of 
the end of tre beam during tne event. 





3 
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By examining these relations, combined with a post-test Inspection 
of the specimens and examination of the films, much about the beam 
response and the failure modes can be detenniricd. The next section 
discusses the features which were common to all tests. 


3.2.2 Descriptio n of Characteristics Common to All Tests 

Six of the relations just discussed for the | ( 15/0/- 15) 5 1 s are 
plotted for eacn of the other 10 laminates tested in f'g. 3.8 through 
fig. 3.37. Three main features of interest ve»*e common to all the 
tests. They are: (1) All the dynamic tests exhibited a third node 

vibratory response; (2) With the exception of two lairirate types, the 
static and dynamic failure modes of the laminates we r e Identical* and; 
(3) All of the beams exhibited varying degrees of binodular material 
behavior. 

Third Mode Response 

Analysis cf the third mode vibratory response is best done through 
examination of the strain time histories. W 4 th the exception of 
lam-nates 9 and 10, all the laminates showed similar responses, 
qualitatively. All laminates showed the initial reverse value of stra-'n 
associated with the initial "W" shape. T hen a vibratory response was 
superposed on a ir.or.cton i cal ly Increasing strain. The vibratory response 
dampened with time and severely dampened after ply failures. The 

frequency of the vibratory response was determ-ned from the experimental 
data for each laminate type and is presented in table 3.1. The 
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Figure 3.8 Static Load and Strain vs. Displacement for the jo]™ 
Laminates ^ 



46 




Figure 3.9 Oynamlc Load and Strain vs. Displacement for the 10] 30 
Laminates 
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Figure 3.10 Dynamic Disalacement and Strain vs. Time for the [01™ 
Laminates 30 
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Figure 3.11 Static Load and Strain vs. Displacement for the 
[(30/0/-30 ) 5 ! s Laminates 
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Figure 3.15 Dynamic Load and Strain vs. Displacement for the 
l(45/0/-45) 5 l s Laminates 
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Figure 3.17 Static Load and Strain vs. Displacement for the 
[ (60/0/-60) S J S Laminates 
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Figure 3.26 Static Lead and Strain vs. DispTacemerr 
(0g/90g) s laminates 
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Figure 3-30 


Dynamic Load and Strain vs. Displacement for the 
(90g/0g) $ Laminates 
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Figure 3.33 Dynamic Lead and Strain vs. Displacement for tne 
L{45/-45/0/90' 4 j s Laminates 
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Figure 3.35 


Static Load and Strain vs. Displacement for the 
[ (G/45/C-45) 3 / 9 O/O/O jy 2 ] s Laminates 


SURFACE STRAIN 


73 


1600 

1200 

800 


(( O/^S/0/^5^/90/0/0^, 


g 0 * 

u 

-400. 
-800 
0 020 . 

0.010 . 

0.000 . 

- 0. 010 

- 0. 020 





l( 0/4S/0/-1 SV30 /O /O 



0 <4 8 12 16 

END DISPLACEMENT. IN. 


Figure 3.36 Dynamic Load and Strain vs. Displacement for the 
[ (0/45/0-45) 3 /90/0/0 l/2 ] s laminates 
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Figure 3.37 Dynamic Displacement and Strain vs. Tine for the 
[(0/45/0-45) 3 / 90 /0/0 1 / 2 l s Laminates 
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theoretical third mode natural frequency for a simply supported beam is 
also presented In table 3.1 for comparison. Reference 6 was used to 


compute the theoretical frequency. Table 3.1 demonstrates that the 
experimental and theoretical frequencies compared quite well. It is 


important to note that even though the initial Impact velocities were 
only 19.6 ft/sec., the maximum strain rates were on the order of 5 
in. /in. /sec. because of the third mode vibratory response. The maximum 


rate occurred during initial snap-through of the third mode "rf" shape. 
If there had been no vibratory response* the maximum strain rate wciild 
have been less. Notice that the strain-time response of laminates 5 and 
7 exhibit vibratory frequencies higher than the third mode, in addition 
to the third mode. 

Laminate 9 failed very soon after Impact in the dynamic tests and 
no sustained vibratory motion occurred. The failure mode was unusual 


and will be discussed in the next section, however, fig. 3.28 shows the 
characteristic Initial reverse value of strain indicating that the 
vibratory response was excited. Due to early failure, *t never fully 
developed. Laminate 10, figs. 3.3C and 3.31, actually did shew the 3rd 


mode response. However, the surface plies failed soon a f ter Impact so 


the response was not measured by the strain gauges. The vibratory 
response is exhibited in the load-time history shown in fig. E.4 in 


Appendix E and it is Quantified in table 3.1. 
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TABLE 3.1 

COMPARISON OF MEASURED VIBRATORY FREQUENCY 
AMO THEORETICAL THIRD NATURAL MOOE 




THEORETICAL 

laminate 

MEASURED 

NATURAL 

NUMBER 

FREQUENCY 

FREQUENCY 


(HZ) 

(HZ) 

1 

227 

237 

2 

223 

231 

3 

205 

207 

4 

176 

178 

5 

155 

158 

6 



7 

152 

149 

B 

203 

206 

9 

* 

254 

10 

112 

120 

11 

172 

177 

12 

187 

185 


♦no oscillatory motion observed. 


Immediate failure after impact 
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Failure Modes 

From examination of the static and dynamic Icad-disulaceirert 
relations, and from the post-test examination of the specimens, it can 
be said that with the exception of two laminate types, the failure modes 
of the laminates under static cr dynamic leading were identical. As had 
beer indicated, the two exceptions to this were the (0g/90g) s and 
(90g/0gj laminates, laminates 9 and 10. These laminates exhibited a 
different dynamic failure mode than was observed with the static 
loadings. A detailed description of the failure modes of eacn laminate 
will be presented in the next section. 


Strain Response 

Examination of the strain response of the laminates revealed that 
at a given load level, a given displacement level, or at a given time, 
the compressive surface strain was always greater In magnitude than the 
tension surface strain. Analysis shows that the compression st“ 3 ir, dee 
to the axial load should have been several orders of magnitudes less 
than the bending strains and effectively tne beam was in a state of pure 
bending. With pure bending the magnitude of the tensile and compressive 
strains should have been the same. For example, in the ((90/0/-90)g) s 
beam, the compressive surface strain was only 5 % greater than the 
tension strain, as shown in Hg. 3.2C and fig. 3.21. On the ether hand, 
fig. 3.11 and fig. 3.12 show that for the ! (30/Q/-30)g] s laminate, 
laminate 3, the compressive surface strain was more than 4G5£ greater 
than the tension strain. Variation of the tensile and compressive 
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strains for the other laminates was somewhere between the 5* and the 40X 
levels. Understanding this phenomena will be important to deve oping a 
failure criteria. Additional tests were carried out to further 
understand this behavior of the beam specimens and the results will be 
discussed in the next chapter. 


3.2.3 Failure Mode Description 

By examining the static and dynamic load-displacement and load 
strain relations, along with observing the movies and the failed 
specimens, the failure modes of each laminate can be described. In this 
section the static and dynamic failure modes will be assumed to be the 
same if no differences are indicated. Figure 3.38 shows a drawing of 
the coordinate system used in the description of the cracking, 
del ami nation, and failure of the laminates. 


Laminate 1 l 0 } 

The unidirectional laninace exhibited a single devastating failure 
event. Figure 3.39 shows a frame of the failure event under dynamic 
loading. Mote the simultaneous failures on the tension and compression 
sides the beam. A similar failure node was observed m the static 
tests. This was the only laminate for which compression fa' ures were 
observed. On the tension side a crack propagated perpendicular to the 
fibers (the y direction in fig. 3.38), across the entire fcidth of the 
laminate, completely failing ten plies. On the compression side a crack 
propagated only partial y across the wdth, in the y direct'on, and 10 
plies deep. The center plies cf the iammate, near the neotra surface. 





Figure 3.39 Fd^ure Mode of Unidirectional Specimen 
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had no fiber failures. In the portion of the laninate which had fiber 
failures, the fibers broomed and splintered. In the portion of the 


laminate with no fiber failures, matrix cracks in the x-z plane 
propagated a^ng the length, presumably initiating on the compression 
side of the laminate due to widthwise Poisson tensile effects. 


Laminate 2 ( (15/0/- 15) 5 ! s 

The load-displacement and the strain-displacement relations show 
that | ( 15/0/— 15) 5 1 5 laminates exhibited from 3 to 5 failure everts. A 
typical failure event consisted of the simultaneous failure of a group 
of 2 to 6 plies on the tension side of the beam. This group of failed 
plies would then delaiminate at least 8 in. along the length cf the 
beam. By examining the failed specimens. It appeared that a crack 
initiated parallel to the fibers in the outer 15 s ply. The crack 
propagated along the 15 1 di rection from one edge of the bean tc the 
other, and it propagated downward in the z direction, to a depth of 1 to 
5 plies, before arresting. As the crack propagated downward, the *ibers 
in the 0 : plies and the -15° plies fractured. So, the 9 C ply and the 
-15° ply failed on a 15’ angle. Exactly now the othe** failure events 
initiated is not clear. However, it is felt they were similar to the 
first failure event. The rema 1 n 1 rg comp r es s ion si de (0* SC n . , * —0 0 
the laminate was undamaged. 


Laminate 3 j (30/0/-3G) 5 i s 

The l (30/0/-30) 5 1 s laminates exhibited 2 or 3 failure events, each 
event involving typical 'y 5 to 10 plies. These ply groups then 
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delaminated at least LO in. along the length of the t>eam. Again it 
appeared that a crack Initiated parallel to the fibers in the cuter 30 r 
ply. The crack propagated along the 30 ; direction from one edge o f the 
beam to the other. In some of the SDecimens the crack turned to 
propagate along the 0 : direction (x direction in fig. 3.3B) for a short 
d-' stance and tnen turret bcc< to the 39" direction. The creek 
propagated dow^a^c in -.no z directicr from 3 to 5 plies, fractur-ng the 
fibers in tie 3" and 30 plies oefere arresting. The remaining 
compression s'cc cf *.h8 . am " rate ! . ’ c 'n. -45X) wjs jrvdiTi^QS d . 


Laminate 4 ! ( 45 >0 - -15 < , ^ • s 


4<though the pcst-fai nscect-c 


* /- • r» r* 4 


♦ hn 

U »C 


i A 

•w VJ J I J| ki" ! > V. 


[ (45/0/-45/)^! s "an '-aces '•evealec nc differences in the failure medes, 
the st r ain responses sho»r in t-g. 3.14 and fig. 3.15 ind'cated 
different oehdv or w € ^ u< -u *j ^ j^cs. ***w u. • -v 

failure inspection of one stat'c and dynamic specimens showed many 
matrix cracks in both the 45 ; and the -45 'amina wn^cn sa^rcurc tre 0 
lamina. F a 1 'a , ‘€ o f tne C : 2 *i proDuT y nit ut?d ■ *?i ^ T5*»« *x w» 
in a 45’ lamina, but the crack direction in the G r ‘'amina followed no 
preferred direction. In some o f the 0 laminae tne crack alternately 
followed 45 c , -45 s and O' directions. In ethers the c r 2ck fcl lowed a 
random jagged line across the width. Examination of fig. 3.16 ircicates 
a unique dynamic strain response for this laminate. The unique response 
is especially evident in the response of specimens 4.4 ind 4.6, 2 nd to 2 
lesser degree in specimen 4.5. Notice that at the point of first ply 
failure, when the tensile strain response terminated, there is no sharp 
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drop In the conp r e$s' ve surface strain, as was seen *n the other 
laminates. Instead, the CQnprp$siyp surface strain remained nearly 
constant for the rest cf the event. ~h ; s suggests that the failure 
necnar^sn *c r th"'s lani r ate could o? different f^om the ones described 
to this point. With this laminate, at the initiation r.f tensile sire 
matrix cracks, which destroyed the tensile strain guaqe, the laminate 
did not fail and inred’ate "y lose load-carrying laminae. Instead the 
tensile side laminae yielded but rpma-npr 'ntact arc continued to carry 
some load as the laminae fai'ed progressively throughout the remainder 
of the event. This type of failure mode for a 45’ angle-ply laminate 
has been reported in the literature (ref. 7), Tb's response was not 
evident in the static data. However, in the dynamic case, the yielding- 
progressive failure response lasted a mere 0.3 e sec, from tne initial 
tensile side matrix cracks, until the end nf the dynamic evpnt: Hyring 
this time interval, the end of the beam displaced 8 in. It is possible 
that in the displacement controlled static tests, the fa"' lure event 
lasted for a similar finite time pe-lod, ever though *_ ho hea ,T wcjid hayp 
displaced less tnan C.C1 ; r . at the quasi-stat'c ciso lacenent rate. So, 
the mode of failure in the static ard cynam-'c tests cod'd have beer tre 
same, even though the ’•eccried strain ’•esponse *as quite di-^p'prt, Pn 
identical static and dynamic failure node Is consistent w'th the 
information learned frcm che post-fa! lure SxdfVnction. 

Laminate 5 \ (60/C/-6C)c ! s 

The [ ( 60 /Q / - 63 ) 5 J c ’d.mirate exhibited d r • r ! 1 1 a major fa! 'ere 

event and a single miner subsequent: eve r t , From examination of *he 
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strain responses, the failure events appear to be instantaneous » !n the 
single major event, an outer group of plies failed. There were many 
parallel matrix cracks occurring in the 60' or -60° plies, and fiber 
failures on the 0 ° plies and souse of the 60 pises. In the subsequent 
minor event the plies failed only partially. This failure mode is 
illustrated in fig. 3.40. A matrix crack occurred in the 60° direction 
and then the crack propagated downward only through the O’ lamina. The 
-60’ lamina remained intact. The -60° fibers required the crack in the 
failed plies to open in extensional and shear modes. The -60’ lamina 
delaminated near the crack but remained attached to the failed plies, 
bridging the crack. 


Laminate 7 [ (90/0/90 ) 5 i s 

The l ( 90 / 0 / 90 ) 5)5 laminate exhibited an initial major failure event 
and a subsequent single minor event. In the major event, the cracks 
Initiated parallel to the fibers in the outer SO' ply. Then the crack 
propagated downward in the 2 direction, fracturing fibers in the 0 C 
plies. The failed ply groups then delaminated 5 to 6 in. along the x 
direction. 


Laminate 8 l(0/90)g| s 

The f(0/90)g|j laminate had only a single failure event. The crack 
propogated across the plies ard downward, similar to laminate 7. One of 
the static and one of the dynamic specimens had some of the 90* plies 
near the neutral axis damaged vrth matrix cracks. There the adjacent C° 
plies remained undamaged. 
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Shear and 
Extentionai Crack 
Opening in 60® 
Lamina 


Intact -60° 
Fibers Bridging 
Cracks 


figure 3.40 Sketch of Partial Failure Mode *n | ( 60/0/-60)g J s 
Lanl nates 
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Laminate 9 {0g/90gl s 

Laminate 9 had a completely different failure mude in the dynamic 
case than in the static case. Statically, specimen 9.1 and 9.2 
exhibited fiber failures, anc matrix cracks in the x-z plane, in the 
outer 0 : lamina groups on the tension and compression sides of the 
beam. Even though this was a static test, as soon as the 0° fibers 
failed, the 90 s core lamina immediately broke apart in chunks and 
scattered about the room. Specimen 9.3, which fail'd at an axial 
displacement level almost 4 ir. less than the others, exhibited no fiber 
failures. However, there were matrix cracks, in the x-z plane, in the 
outer 0* lamina groups and the 90 r lamina again dispersed about the 


room. 

In contrast, under dynamic loads there was no failure in the outer 
(T lamina groups. Soon after the impact event, a del ami nation initiated 
between the group of O’ lamina on the compression side of the bean and 
the 90° core laminae. The delamination propagated along the entire 
length of the bean. The 90 l core laminae remained attached to the C“ 
lamina on the tension side of the beam. "The laminate acted as two 
separate thin beams, figure 3.41 shows a sequence of four frames at the 
initiation of the del ami rat ion. Notice that the delaminatior initiated 
between the center peak of inital “W" shape and the inflection point of 
that peak, i.e., about at the quarter span near the top of the beam. 
The construction of laminate 9 produces Interlaminar tensile c 2 stresses 
at the free edge when the 1 aminate is placec under tension. These are 
caused by the Poisson ratio mismatch of the 90' core and the 0° outer 
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Figure 3.41 Dynaaic Failure Mode of (0g/90g) s 
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lamina group. The top half of the plies In the laminate were under 

tension in the center peak region. In addition, due to the deformed 

shape, the t x2 stresses reach a maximum at the inflection point of the 

center peak. And due to the high amplitude lead spike shown In fig. 

3.27 the stresses are particularly high at this time in the event. 

Although the stresses were similar - for the other laminate types, 

laminate 9 had clustered plies. The interlaminar z 2 stresses are much 

greater in orthotropic laminates which contain clustered plies. Only in 

the dynamic testing of laminate 9 did large values of both 

the t ard c stresses occur simultaneously. Therefore It was thought 

xz z 

that the combination of the o z and - %1 stresses initiated the 
delamination in this particular laminate type. 


Laminate 10 I90g/0gl s 

Lamlrate 10 exhibited differences in the failure modes between the 
dynamic and static loadings. Under static loading matrix cracks 
appeared In the 90° plys on the tension side with spacing of 1 to 2 
inches. However, the 90" lanrnae remained attached to the 0 e core. Ho 
matrix cracks appeared on the compression side. Thus the 3C' laninae cn 
the compression side contributed to tne laminate stiffness. Notice the 
unusually hi^h strain values of over 2% in fig. 3.29. This exceeded the 
range of the recording apparatus. However under dynamic leading doth of 
the 90' outer lamina groups separated in chunks from the O’ core. 
Therefore the bean acted as a single thin undirectional beam. This 90’ 
lamina separation was probably due to the inertial force and reverse 
curvature from the excited vibratory motion. 
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Laminate 11 I (45/-45/0/90) 4 l s 

The strain-displacement relations for laminate 11 showed that It 
exhibited behavior similar to laminate 4. However with this laminate 
the failure mode was not as distinctive. Some of 'he dynamic specimens 
had several minor sharp drops in the compressive strain value. In 
addition to the plastic hinge behavior. 

Laminate 1? I (O/45/O/-45) 3 /9O/O/0 1/2 ! s 

Laminate 12 failed with one failure event, exhibiting Instantaneous 
ply failures and longitudinal del ami nations. The ply failures 
propagated more or less in a perpendicular fashion across the beam, but 
with no preferred direction, similar to a paper tear. Figure 3.42 shows 
a frame from the failure event of this laminate type under dynamic 
loading. 


3.2.4 Comparison of Static and Dynamic Tests 

One of the most dramatic comparisons between ;he static and dynamic 
results concerns the end displacement at initial failure. For virtually 
all laminates, the value of end displacement at failure was less for the 
dynamic cases than it was for the static cases. Likewise, the strain 
level at failure was always lower in the dynamic cases than in the 
static cases. The strains are, of course, related to the end 
displacements. Even though inertia and toe third mode oscillation make 
the relation less direct for the dynamic case than for the static case, 
it is felt to be significant that strain rate or displacement rate 
appears to influence failure strains. 
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Figure 3.42 Failure Hode of ! (0/45/0/-45) 3 /90/0/ 1/2 ) s 
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To further compare the static and dynastic events , their end load 
vs. end displacement relations were Integrated. These Integrated 
relations are presented in fig. 3.43 through fig. 3.48. These Integrals 
are associated with work done by the force at the top of the Ossiu. This 
type of relation was explored as a method of determining if there were 
any energy absorbing mechanisms displayed In the dynamic response that 
were not evident in the static response. In addition, the integral 
effectively smoothed the dynamic load repsonse which allowed for 
comparisons to the already smooth static load-response. 

On the figures the vertical axis gives the *or< required to 
displace the beam an amount given on the horizontal axis. The points on 
the work -displacement relation where the slopes abruptly change 
correspond to displacement values at which failure occurred. Again 
laminates 9 and 10 showed unusual differences in the static and dynamic 
responses and will be discussed separately. 

As a group, all laminates except lasj’nates 9 and 10 showed similar 
static and dynamic end load work vs. displacement relations. However, a 
close examination of the relations reveals there were some differences 
between the static and dynamic responses. This can be seen by examining 
fig. 3.43, for example. For the first inch of end displacement, the 
static and dynamic relations are divergent. The work required to 
displace the beam axially through the first 3/4 in. is greater in the 
dynamic tests than in the static tests as a result of the inertia of the 
beam resisting deformation more in the dynamic case. However, after the 
motion begins, it appears that overall, the work required to displace 
the beam Gynarr'ca 1 ly a given -'ncrement is equivalent to the work 






END WORK, 1N-LBF ENO WORK. IN-LBF 
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Figure 3.44 End Load Work vs* Displacement for the ( (30/0/-30)gJ s 
and { (45/07-45)5 ] $ Laminates 
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Figure 3.45 End Load Work vs. Displacement for the [ (6O/0/-6G) 5 ] 
and l (90/0/90)gl s Laminates 
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Figure 3.47 End Load Work vs. Displacewnt for the i\90g/0 8 0] $ and 
{ (45/-45/C/90) 4 I s La*inetes 
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required to displace the beam statically ths sams increment. Tf» 1 s 
indicates that the displacement rate of doing work is comparable for 
Doth the static and dynamic everts. The other laminates in this group 
show comparable behavior up to first failure. The work required to 
displace the specimens dynamically 3/4 in. was fairly consistent ^rom 
specimen to specimen and was about 400 in. -lbs. The work required to 
statically displace the beams 3/4 in. varied with beam stiffness. 
Therefore, the softer laminates, such as 5 and 7, showed greater 
divergence between the static and dynamic events, as can be seen in fig. 


3.45. 

An interesting general conclusion concerning work can be made by 
further comparing the work for the static and dynamic cases. As stated 
earlier, failure of the dynamic specimens occurred at a lower value of 
end displacement than failure of the static specimens. As a result, the 
static specimens required more work to displace the full 16 in. than the 
dynamic specimens did, even though the dynamic specimens required 
significantly more work to displace the first 3/4 in. 

Finally, fig. 3.46 dramatically shows the differences in the static 
and dynamic failure modes for the [0g/90gj laminates, laminate 9. The 
initial dynamic response of laminate 9 is compa r ab1e to the other 
laminates, i.e., requiring 400 in. -lb. of work to move the end 3/4 In. 
However, soon after the impact event the beams delaminated and separated 
into two thin beams. The two thin beams required considerably less wo r< 
for deformation. There were moderate differences Detween the static and 
dynamic failure modes for laminate iO, the (9Gg/0gi s laminate. Srortiy 
after impact, the outer plies separated from the beam and t h e beam acted 
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as a single thin unidirectional beam. In the static test, the outer 
plies cracked but remained intact and contributed to the stiffness of 
the laminate. Therefore, the slope of the static work-deflection 
relation was different for the static case than for the dynamic case. 

3.3 Suwary of Experimental Observations 

Laminates 9 and 10 exhibited unusual and distinct failure modes 
under both static and dynamic loadings. However, because the laminae in 
these beams were clustered, these laminates are not of practical 
interest. In fact they were designed to nurture differences in static 
and dynamic loadings. In this they were successful* The rest of the 
laminates exhibited nc differences between the failure modes under 
static and dynamic loading. 

The progression of the failure* or the number of failure events, 
varied from six events, for specimen 2.3, to one event. Also, in 
general, the specimens that were tested dynamically fai 'ed at strain 
values and axial displacement values less than those of the static 
specimens. These results will be explored further in ch. 5. 

Finally, variation between the tensile and compressive surface 
strains were observed in the static and dynamic tests. The variations 
ranged from an increase in compressive strain of from St above the 
tensile strain to increases of up to 4 0%. These phenomena will fee 
investigated and discussed further in the next chapter. 


Chapter 4 

INVESTIGATION OF DIFFERENCE BETWEEN COMPRESSIVE 
AMD TENSILE MODULI 


Examination of the static and dynamic strain response of the 12 
laminate types showed that the compressive side surface strain at the 
center of the beam was greater than the tensile side strain. This 
difference in tie surface strains varied from S% for laminate 7 to 40* 
for laminate 4. At first It was felt that the greater compression 
strain was due to there being a net compressive axial force on the 
laminate. However, calculations show that these compressive strains can 
be effectively ignored since they are always at least 2 orders of 
magnitude less than the surface strains actually measured. In addition* 
there is no shear deformation at the center of the beam due to the 
symmetry cf the loading. So, the center section of the beam is 
effectively in a state of pure bending. With pye bending and a linear 
elastic material, the magnitude of the tensile and compressive strains 
should be identical but this was not the case. The observed behavior 
was known to be elastic. Recall, It was found in the static tests that 
the loading and unloading strain response curves were coincident, unless 
there were ply failures. Therefore the difference in compressive and 
tensile strains was due to eithe* - a nonlinear elastic effect or the 
laminates were perhaps exhibiting b^modular elastic behavior. 

To determine which of these two phenomena were occurring, laminates 
l through 7 and laminate II were se'ected for further study. Further 
tests were conducted to empirically determine the effective compressive 
and tensile bending modjli. 
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4.1 Experimental Set Up 

Since the unusual material behavior was observed under static and 
dynamic loadings, for simplicity, these experiments were conducted in a 
static test apparatus. To measure the elastic properties. In addition 
to the strain data, knowledge of the Internal bending moment is 
required. The bending moment at the center of the end loaded beam Is 
equivalent to the end load times the distance between the line of action 
of the end load and the neutral surface of the bean at its point of 
maximum lateral deflection. The maximum deflection occurred at the 
center of the beam. This is illustrated In fig. 4.1, D being the 
distance of Interest. To fad late measurement of the lateral 
deflection, a modification to the original static test apparatus was 
made. A rigid bar was mounted perpendicular to the moving crosshead and 
parallel to the undefomed beam. Tnls bar provided a reference from 
which the lateral deflections could fee measured . A1 1 other aspects of 
the set Lp were identical to the set uo for the static tests performed 
earlier. However, different strain guage amplifiers and a different 
load frame were used. 


4.2 Data Aquisition 

As in the previous static tests the load was ne asured by the 1000 
pound load cel’ and the end load-end displacement relation was recorded 
on a chart recorder. The strain guage signals were conditioned with 
signal amplifiers. Strain-load relations were recorded using an X-Y 
plotter. In addition, the axial end load, axial end displacement, 
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center lateral displacement, and the two strain values were manually 
recorded from digital voltmeters. 

4.3 Test Procedure 

first the beam specimen was placed in the load frame and the strain 
guages were connected. Then the load cell was zeroed and the strain 
guage amplifiers were balanced. The beam specimen was displaced axially 
and monotonically in 1/2 In. increments up to 85* of the axial 
displacement value at which failure was first seen in the previous 
static tests. At each 1/2 in. increment the lateral deflection was 
measured by hand with a ruler, then recorded. At this time the axial 
displacement, the end load and the two strain values were also 
recorded. This procedure was repeated until the predetermined axial 
displacement was reached. Then the displacement direction was reversed 
and the beam was unloaded until the beam reached the zero displacement 
point. The strain- load relation for unloading was checked to see if it 
was coincident with the loading relation. 


4.4 Data Reduction 

The data was reduced to determine the tensile and compressive 
moduli of the laminate. The reduction procedure utilized the fact that 
the beam was in static egui librium. In addition, three assumptions were 
made. The first assumption was that the material properties could be 
smeared through the thickness. In other words, the stress was assumed 
to vary linearly from the outer surface to the neutral surface. The 
second assumption was that there was a linear variation of strain 
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through the thickness of the bean, i.e., the Klrchhoff assumption was 
valid. This assumption should be valid In that the bean was effectively 
in pure bending. Finally, the third dSSuiMptlon WS that the T*s«1 fiats 
was bilinear elastic. This may or may not have been the case. If It 
were the case then the determined moduli should not vary wltn Increasing 
curvature, load, or strain. If the moduli did Indeed vary with 
curvature, load, or sttain, then it could be concluded that the laminate 
was not bimodular but rather exhibited some other form of nonlinear 
elastic behavior. However, the precise nonlinear relation could not be 
determined from the analysis. 

From the Kirchhoff assumption, knowledge of the total beam 
thickness and the empirical values of surface strain, the location of 
the neutral surface within the crosssection could be found. The 
location was given by 


t 


1 



and 



(4.1) 


(4.2) 


where 

tj * distance from the neutral surface to the tension surface, 
t£ = distance from the neutral surface to the compression surface, 
T - total beam thickness, (mid thickness. Table 2.2) 
z c = strain at compressive surface, (absolute value) 

- strain at tensile surface. 
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Figure 4.1 illustrates the geometric quantities involved. By 
assuming the net axial stress resultant could be set to zero, and by 
using the bldiaterial assumption, 3 relation between the tensile and 
comresslve moduli could be found. That relation is 

| E t e t t 1 b = | E c e c t 2 b , (4.3) 

where 

E* - unknown tensile modulus 
€ c = unknown compressive modulus 
b * beam width (* 2.00 in.). 

Finally by enforcing moment PCjw ^11 br IImm for the upper or lower hal • of 
the beam. It is fojnd, 

PD = E t e t t 1 b) (§ tj} * \ £ C c°t 2 b) (§ t 2 ) (4.4) 


where 

P = end load 
D - moment arm 

The values of ti and t 2 can oe determined from eqn. 4.1 and eqn. 4.2. 
Solving eqn. 4.3 and eqn. a. 4 for the moduli yields: 


E c . 3P0 

e c (t 2 + tjt 2 'b 

F t , 3P0 

i t ? * 

e(tf + V 2 > 


(4.5) 


(4.6) 
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^ c t 

Multiplying both sides of egn. 4.5 find sqn. 4.6 by e find c , 

respectively, yields equations for the surface stresses, namely, 


3 - 


3P0 

;t 2 + t l t 2^ 
3P0 

^ + W 


The beam curvature can be computed from 


< = 


*1 * 


(4.7) 

(4-8) 


(4.9) 


Using this procedure the compressive moduljs, the tensile modulus, the 
stress at the tensile surface, the stress at the compressive surface, 
and the beam curvature were computed at each displacement increment. 
The raw data used for the computations are presented in Append* x 6. 


4. 5 Results 

The modul i -curvature relaticrs are potted *or each laminate tested 
and are given in fig. 4.2 through fig. 4.5. The horizontal axis shows 
the beam curvature and the vertical axis shows the emp^nically 
determined bending tensile arc compressive moduli. All the laminates 
exhibited an elastic response. This was known because the graphically 
recorded load-strain relations wene coincident for loading and 
unloading. The uni directional laminate, fig. 4.2. shewed the 
compressive modulus to be 235f less than the tensile modulus at high 
curvature levels. Published values for the fiber direction modulus 
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indicate that compressive modulus should be only 9.75X less than the 
tensile modulus. As shown in fig. 4.4 and fig. 4.5, laminates 6 and 7 
exhibited the least amount of ulmodularity. Also, laminates 6 and 7 
exhibited nearly constant values of the modulus with curvature, or more 
specifically, a linear response of stress with strain. In the other 
laminates the compressive modulus decreased with curvature and the 
tensile modulus increased with curvature. This is especially evident in 
laminates 4 and 11 shown in fig. 4.3, and fig. 4.5, respectively. 
Because the moduli were not constant with curvature there appears to be 
a nonlinear elastic stress strain response more complex than 
bimodular. Laminates 3 and 4, shown in fig. 4.3, exhibited a remarkable 
amount of bimodulanty. Notice that for both of these lami nates the 
compressive modulus was nearly 50* less than the tensile mcculus at high 
curvature levels. 

Effective stress-strain re'ations fo r the different laminates were 
also determined from the bending tests. Though it is known that the 
stresses change value dramatically from one lamina to the next (in 
contradict ior to one of the assumptions used to reduce t “e data), a 
stress-strain relation for the laminate-as-a-whole in tension and a 
stress-strain relation for tne laminate-as-a-whole in compression we^e 
computed. 8y us^ng the measured values cf 2^ and t and the stresses 
from eqn. 4.7 and eqr. 4.8, effective stress-st r ain relations we^e 


empirically determined. Figures 4.6 through 4.9 show these stress- 
strain relations. It is luportar t to review the derivation cf eouat 1 cns 


in this chapter to put the stress-strain benavior snown In fig. 4.6 
through fig. 4.9 ^nto context. However, it appeals that, to a f'rst 
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approximation, the laminates did exhibit bimodular behavior. A close 
examination of the figures reveals that the two portions of the stress- 
strain relation are not perfect straight lines. Therefore, a ircre 
complex nonlinear elastic behavior may be the case. 


4.6 Error Analysis 

The compressive and tensile moduli were determined from 5 measured 
quantities. From eqn. 4.5 and eqn. 4.6 

E c = E c (P,0,c t ,c c »TJ (4.10) 

and 

E l = E t i'P,0,g t ,£ C ,T j (4.11) 

Uncertainties in each one of the five independent variables contribute 
to the overall uncertainty in the moduli. The uncertainty in t*e *oad 
was 1% of the maximum value of the load cell, in this case 10 lbs. The 
lateral deflection, which was measured by hand, had an uncertainty of 
0.05 in. There were many factors contributing to the uncertainty in the 
strain measurements. They were: (1) The finite distance between the 

foil grid of the strain guage and the surface of the beams (2) Drift in 
the amplifiers, (3) The transverse sensitivity of the strain gauge, and; 
(4) Misalignment of the guages or the specimen. It was felt that all 
these errors together contributed to an uncertainty of .00010C in. /'in., 
or *1*, of the maximum strain measured. Finally the uncertainty in the 
thickness measurement was taken to be 0.0025 in. or 2%. This accounts 
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for the uncertainty In the measurement and the variation of thickness 
across the width of the specimen. The way in which each of these 
uncertainties contribute to the uncertainty in the computed moduli can 


be given by the relation 
, I t jE . ,'li U . (ll 




(4.12) 


where is the uncertainty In the moduli and Wp, ...» W T are the 
uncertainties in the individual variables. 

The uncertainty calculations were performed «ith each calculation 
cf the compressive and tensile moduli. The results of the uncertainty 
analysis shew that the uncertainty In the moduli calculations can be 
quite high, 2-3 MSI for the first data ooints taken at lew values of 
axial displacement and "ad. However, the uncertainty i fi the ImGuu ^ 1 IS 
much lower, 0. 5-1.0 MSI, for the data taken at large values of axial 
displacement load. The uncertainty of the moduli were most sensitive to 
the uncertainty in the strain iReasurericnts. The uncertainty in the 
strain measurements were the most difficult to quantify. If a vertical 
error bar had been pin ,ed with each data point in the moduli -curvature 
relations of fig. 4.2 through fig. 4.5, a horizontal tensile •nodulis vs. 
curvature relation and a different horizontal compressive modulus vs. 
curvature relation could have been drawn within the error bars. This 
would mean that the moduli did not vary with curvature, and hence load 
or stress level. Therefore, the apparent nonlinear behavior other than 
bimodular is not conclusive. However, it appears there is still 
substantial experimental evidence of a difference between the tensile 
and compressive moduli. 
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4.7 Observations 


Reasons for the nonlinear behavior are rot entirely clear. !t is 
surmised that in tension the fibers align themselves with the principle 
stress directions and in compression the fibers buckle, or turn away 
from the principle stress directions. These Uc tendencies result In 
different stress-strain behavior in tension than in compression in a 
bending situation. This behavior was entirely unexpected and the 
thorough examination of it was not with in ths overall $C3l s of this 
study. However, the flexural moduli empirically derived here will be 
used in ch. 6 with a finite element program to predict the deformation 


resporse of the beams. 



Chapter 5 

INVEST IGATIOM OF THE FAILURE RESPONSE 


Predicting the response of an aircraft fuselage to crash conditions 
involves both the prediction of the dynamic deformation response of the 
structure, and the prediction of the on set and extent of failure in the 
structural elements. This chapter focuses on tne later by investigating 
the experimentally recorded ultimate surface strains in each of the 
laminates. Because of the uncertainty of the material behavior, i.e., 
blmodular vs. some other form of nonlinear elastic behavior, the onset 
of failure in the laminates was investigated using empirically measured 
strain values. This is in contrast to using stresses. The static data 
was studied to investigate the dependence of ultimate strain values with 
the laminate type. Then the static data were compared to the dynamic 
data to determine the effect of the dynamic response on the ultimate 
strain of the laminates. Since the failure modes for laminates 9 and 10 
were significantly different in the static and dynamic tests, the 
ultimate failure strains could not provide a meaningful comparison for 
these laminate types. Therefore, tne strengths of laminates 9 and 10 
were rot investigated in this chapter. 

5.1 Oata Reduction 

The ultimate longitudinal tensile and compressive 
strains, e* and were determined f rom the experimental strain- 
displacement response of the static and dynamic tests. Tne value 
of was taxer to be t h e 'ongtitudinal tensile strain ya'ue at the 
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displacement just prior to guage failure. The value of was the 
compressive strain at the same displacement. The ultimate strains for 
all of the static and dynamic specimens tested are presented in Appendix 
f. The average ultimate e* and were computed for the static 
specimens and the dynamic specimens of each laminate type. The strain 
transformations were used to transform the average values of 

longitudinal strain into the material pr >rcple strains ^2 “"** ^12* 
Since al 1 of the failures (except in the unidirectional laminate) 
occurred on the tension side, only the tensile ultimate strains were 
transformed. In addition, the maximuss fiber strain in the C Isnsina 
closest to the tension surface was computed. This ply was of Interest 
because it was the most highly stressed lamina, though not the most 
highly strained. A summary of all the computations IS presented m 
table 5.1. 


5.2 Quantitative Characteriza ti on of Failure 

The mechanics of failure in a composite beam under bending loads is 
complex. For the laminate types selected for this study, failure was 
defined as the initial loss of load carrying capacity. This was usually 
associated with fiber failure in one or more lamina. As described in 
ch. 3, the failure events probably occurred in the fo 1 lowing order. 
First, a matrix crack initiated in an angle ply lanina which was located 
on or near the tension surface of the beam. It is speculated that the 
crack then caused an axial and/or shear strain concentration in the 
fibers in the ply adjacent to the crack. This is depicted in f-g. 
5.1. The initiation of the matrix crack *n the outer most argle ply and 
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TABLE 5.1 
FAILURE STRAINS 


Liainata 

fcjcfcer 

Avg 

#■ 
c ~ 

UH 

(*) 

Avg 

c 

c 

Ult 

(*) 

c. 

a 

U’t 

(*) 

Tension 

*2 

'Jit 

(*) 

*12 

m 

(*) 

T ersor 

Po»y 

In C* 
la*. 
MirtSt 

ten. tu^ace 

(*) 

Strain 

Cone. 

factor 

Static 1 

1.6ls.il 

1.611.11 

1.61 

-.48 

C.00 

1.42 

1.61 

1.00 

Dynamic 1 

1. S0t-04 

1.701.02 

1.50 

-.45 

o.oc 

1.36 

1.5C 

l.OC 

Static 2 

0.9S;. 01 

0.991.04 

0.75 

•.50 

-.72 

0.8$ 

0.79 

2.04 

Dynaalc 2 

0.7S:-00 

0.881.00 

0.66 

-.44 

-.64 

0.66 

0.7C 

2.14 

Static 3 

1.02:. 02 

1.37;. 03 

0.55 

-.55 

-1.81 

1.58 

0.94 

1.71 

Dynaalc 3 

0.93;.01 

1.291.01 

0.45 

-.55 

-1.65 

1.41 

0.86 

1.74 

Static 4 

l.lli.01 

1.521.03 

0.17 

0.17 

-1.61 

1.20 

l.o: 

1.59 

Dynaalc 4 

l.OOt.Ol 

1.401.01 

0.15 

0.15 

-1.69 

0.96 

0.92 

1.63 

Static 5 

1.181.06 

1.431.06 

0.02 

0.79 

-1.33 

2.93 

1.09 

1.46 

Dyn»lc 5 

0.991.03 

1.191.04 

0.02 

0.67 

-1.12 

2.17 

0.93 

1.61 

Static 7 

1.331.01 

1.4s, 01 

0.05 

1.33 

O.OC 

6.47 

1.2* 

1.30 

Dynaalc 7 

1.221.03 

1.331.03 

0.04 

1.22 

o.oc 

5.51 

1.14 

1.32 

Static 8 

1.57i.04 

1.741.04 

1.57 

3.07 

o.x 

1.01 

1.57 

1.02 

Dynaalc 8 

1.441.05 

1.601.07 

1.44 

0.06 

o.x 

0.98 

1.4* 

1.04 

Static 11 

1.381.00 

1.701.00 

0.48 

0.48 

-1.80 

1.96 

1.19 

1.35 

Dynaalc 11 

1.211.02 

1.511.02 

0.42 

0.42 

-1.58 

1.56 

1.03 

1.46 

Static 12 

1.341.37 

1.661.09 

1.34 

0.24 

O.X 

1.04 

1.34 

1.20 

Dynaalc 12 

1.291.06 

1.671.10 

1.29 

0.23 

O.X 

1.07 

1.29 

1.16 
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Figure 5.1 Strain Concentrators Near Matrix Ci*a Cfc 
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the failure of the fibers in the adjacent ply fissy or «5sy not nave 
occurred simultaneously. 

To examine failure In the laminate, a phenomenological approach and 
a mechanistic approach were taken. The phcncricrc leg i CAl Approach tfiksn 
was the strain tensor polynomial (ref. 8). This approach treats the 
composite as an anisotropic but homogenoeus material. The strain tensor 
polynomial interprets failure as the occurence of any definable 
discontinuity in the material response. This interpretation of failure 
may not correspond to loss of load carring capability. The mechanistic 
approach taken was to examine the strain in the 0’ loini na on the tension 
surface, or the 0" lamina nearest to the tension surface. Because this 
was the most highly stressed lamina, failure of this lamina would 
initiate the loss of lead carrying capability. 


5.2.1 Strain Tensor Polynoroinal 

The simplified plane-stress version of the strain tensor 


polynomial can be expressed as 



4 

where E^ and E^ a^e the absolute value of the ultimate extension and 
compressive norma 1 strains in the 1 direction, i = 1, 2, end G is the 
ultimate shear strain. The interaction term was taken to be zero. The 

i 

values of E^ and E^ used for the static specimens were the ultimate 
static strain values of the unidi^ect^cra 1 laminate, laminate 1. 

i 

Likewise, the values of Ei and used for the dynamic specimens were 
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the ultimate dynamic strain values of the dynamics m> tested 
unidirectional lam-hate. The other values used were: » 

0 . 47 *, = 0 . 67 * and S = 1.93*. These were taken from published 

results *or AS4-3502 and were used fcr both the stat : c and dynamic 
spec hers. The ultimate strain values were used in the tensor 
polynomial to predict failure. If the value of tne tersor polynomial 
was less than 1 at failure, the polynomial -oj'd have been 
nonccnservative in the prediction of failure. If the value of the 
polynomial was greater than 1, it would have been conservative In the 
prediction of the loss of load carrying capability. However, in the 
later case the polynomial may have accurately predicted the occurence cf 
material failures, such as matrix cracks in the outer layer of the 


laminate, even though it conservatively predicted tne occcrance of the 
loss of load capacity. T ne results of the tenser polynomial 
calcinations are presented in table 5.1. 


5.2.2 Mecha n istic Approach 

For each laminate type, the maximum fiber strain ir the 0 lamina 
nearest to the tension surface was calculated and is presented in table 
5 2 results Show thdt this f ; £>€ ** St r 2H n wd$ n Ot CO n Std^t 
with laminate type. Thus this simple maximum fiber strain criteria 
could not oe used as an indicator cf laminate behavior for more general 
laminates. However, the f^ber failures it these 0" lamina actually 
occurred in a region of strain concentration, under a matri< crac< in an 
adjacent lamina. This was shown in fig 5.1. ’he ultimate tensile 
strain in the finer direction for tne material w a$ known frpm the 
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statically and dynamically tested undirectlonal laminate. Therefore, 
for each laminate type, the strain concentration factor associated with 
the crack in the adjacent lamina could tse determined. This sias dens by 
dividing the laminate ultimate strains (again as average str; in) by the 
ultimate tensile strain of the material in the fiber direction. T he 
strain concentrations computed this way are presente 


*l I (U J* l* • lie 


strain concentration factors for the static tests used the tensile 
strain at failure in the static tests of the unidirectional specimens. 


i.e., 1.61*. The strain concentration factors for the dynamic tests 
used the tensile strain at failure in the dynamic tests of the 
unidirectional specimens, i.e., 1.50*. T nus, for example, for laminate 
4, the static strain concentration factor that the cuter 45’ ply caused 
the adjacent 0 : ply to experience was 1.612/1.01% = 1.59. Similarly, 
for the dynamic case, the strain concentration factor was i.5G%/'0.92* * 


1.63. 


5.2.3 Results 

Upon examination of the results it should be noted that there '$ c 
consistent difference between the static and dynamic resjlts. The 
dynamic failure strains are d : wdy$ 'ower tnan tne static fai'ure 
strains. Except for laminate 12, the computed strain concentrations are 
always higher and the tensor polynomial value is always lower for the 
dynamic case. Reasons for this are discussed in tne next section. 

Examining the results of both the tensor polynomial and the strain 
concentration factor provides a qualitative understanding of the effect 
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of laminae orientation on the ultimate strain in the : animate. for 
the l(0/O/-e)J family of laminates, laminates 2 through 7, both the 
static and dynamic 0 3 fiber concentration factors decrease with 
increasing off-axis angle. Also the tensor polynomial value increases 
with increasing off-ax^s angle. for example, in the static case of 
laminate 2 which has a 15% off axis fiber angle, the tensor polynomial 
had a value of 0.85 and the strain concentration factor was 2.04. The 
low value of the tensor polynomial could be contributed to the 
generation of interlaminar shear stresses near the free edge of the 
laminate which were not accounted for in the simplified plane-stress 


version of the tensor polynomial. The shear stress is especially 
high in laminate 2 due to the large n x mlssmatch between the 15 c and 
C° laminae, and the 15° and -15° laminae. The low value of the tensor 
polynomial suggests that the initial material discontinuity in the 
laminate ana the loss of load carrying capacity of the laminate occurred 
simultaneously. This is consistent with the high value of strain 
concentration calculated for this laminate. The initial 15 matrix 
crack in the surface ply of the laminate might have caused high shear 
strains as well as nigh tension strains at the base of the crac* «.he 
adjacent O' lamina. There is evidence* however, that the shearing 
effect causes higher strain concentrations than the tension effect. It 


is speculated that the amount of shear strain concentration is related 

to the n tern of the angle ply lamina. For the ( (e/0/-e) family, 

the computed strain concentration factor is a maximum for the 15* angle 

ar.d so is the value of - . The computed strain concentration reaches 

xy 9 x 

winimums at 0 s and 90° as does • . A $ stated above, laminate 7 , with 

xy 9 x 
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off-axis angle plies of 90°, had the lowest value of strain 
concentration of laminates in the {{5/0/-*)g! family* It also had the 
highest value of the tensor polynomial. This suggests that the 90’ 
outer ply cracked much before the load bearing 0’ failed and a SG r crack 
had a less severe effect on the adjacent 0° than did an off-axis crack* 
Comparing laminates 7 and 8 and also comparing laminates 4 and 12 
snows that it is more severe for a O' lamina to nave cracked angle plies 
on both sides of it thar if there is a cracked angle ply to one side 
only. Laminate 7 had 90’ laminae or both sides of the 0* lamina nearest 
to the tension surface. The static strain concentration factor was 
1.30. Laminate 8 had a 0 s lamina on the surface with just ore adjacent 
90* lamina. There the strain concentration factor was 1.02. Laminate 4 
had 45° laminae on botn sides o' the 0’ lamina nearest to the tension 
surface. The strain concentration factor was 1.59. Laminate 12 had a 
0“ lamina on the surface with just one adjacent 45' lamina. The strain 
concentration was 1.20. These numbers indicate that cracks in the 
laminae on both sides of a 0° lamina are more detremental for the 0* 
lamina than if only one of the adjacent laminae is cracked. 


5.3 Differences in Static and Dynamic Failure 

As stated earlier, the results show that the ultimate strains in 
the dynamically tested specimens were lower than the Ultimate strains of 
the statically tested specimens. Two reasons are speculated. First, 
the composite may have strain-rate dependent material properties. 
Second, the dynamic specimens undergo a more severe lead history. 


t 

* 


maximum strain rate in the dynamic tests was approximately 5 
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in./ln./sec. A reduction ir the ultimate strain of matrix dominated 
failures might be expected at this Strain rate. However, a reduction In 
the ultimate strain was seen for fiber dominated failures as well. It 
is felt that strain-rate effects for fiber dominated failures are not 
likely. Therefore it was felt that other mechanisms must contribute to 
the lower ultimate strains in the dynamic tests. The dynamic specimens 
experienced different leading histories than the static ones 5 Recat 1 
that in the dynamic tests there was an initial high-amplitude load spike 
and an initial "W" deformation shape of the beam. The *W* deformation 
was such that the lamina which failed under tension were initially under 
compression. Also the dynamic specimens were subjected to a third mode 
vibratory response which could have fatigued the specimens. The dynamic 
loadings could have caused damage in the lanlrate on the "icromechanicai 
level early in the loadirg history and these could have contributed to 
the lower ultimate strain levels later in the 'oading history. 

5.4 Prediction of the Extent of Failure 

As a final note, this chapter has investigated the conditions which 
lead to the onset of failure. Failure was defined as the initial loss 
in load carrying capability. However, predicting the amount of drop in 
the load upon failure is as important as the prediction of the onset of 
failure. Recall from the lead displacement relations discussed in eh, 3 
that laminate 2 had from 3 to 6 failure events which resulted in 
relative small drops in load when a group of 2 to 6 plies failed 
simultaneously. Laminate 3 had 2 or 3 failure events, .where 5 to 10 
plies failed simultaneously. All the other laminates had 1 major 
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failure event, where 40X to 60< of the plies in the laminate failed 
simultaneously, producing a drop in the load of up to 80*. Further 
investigation is needed to determine the relation between the extent cf 
failure and laminate type. Such studies are Important if the collapse 
of a composite structure is to be predicted. During collapse, some 
structural members will fail, others wi 11 partially fail, and some will 
remain intact with load being transferred and retransferred throughout 
the structure. Knowing how to predict partial loss of load capacity of 
the individual elements would be important for predicting the collapse 
behavior. 


Chapte** 6 

PREDICTION OF DEFORMATION RESPONSE 

To predict the deformation response of the beam, an existing finite 
element program was used. The static load-deflection response *ss 

predicted for each of the laminates and was compared with the 
experimental results. For laminates 1, d, and 7, the dynamic test 
conditions were modeled. The finite-element predictions for the load 
deflection, load-time, and the displacement-time were compared with the 
experimental results. Because of the expense in running the computer 
program, only 3 of the 12 laminates were selected for the dynamic 
analysis. 

6.1 Finite Element Program 

The finite element program used was capable of computing a 
nonlinear transient response of a structure subjected tc time varying 
loads (ref. 9). The program allowed for large geometry changes by 
using a co-rotational coordinate system in the deformation mcce" . The 
analysis is preformed by direct minimization of the scalar energy 
function. 

6.2 Material Model 

The program was written for crash analysis of aluminum 
structures. Therefore the materia 1 model was designed for isotropic 

linear elasto-plastic materials. Although nonlinear elastic materials 
cannot be modeled, linear bimcdular materials can be modeled and were 
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used for the analysis. In Addition* the program cannot ROdel 
anisotropic materials. However, since the beam is one dimensional and 
the axial loads are insignificant compared to the bending loads, the 
compressive and tensile flexural moduli of the laminate ~2S ussd in the 
material model, which raises another issue. Since the material was 
assumed bimodular, the neutral surface of the beam was not coincident 
with the midplane. The program has the capabl 1 1 tj( to compute the 


position of the neutral surface at every position along the beam and at 
every time step. However to do so Increased the computation time by a 
factor of ten. So this option was net evoked for the analysis performed 
in this chapter. However, the effect of tracking the location of the 
neutral surface on the static response is examined In Appendix H. As 
shown in the appendix, the error 1$ significant, especially for laisinats 
3. However, for the dynamic case one computer run using the neutral 
surface computations would cost SiOOO. In addition, as will be seen 
shortly, there were other problems with the computer results. Thus the 
error was tolerated. 


6.3 Material Properties Used 

To perform a successful analysis, accurate material properties are 
required. Determining the properties of a laminated composite material 
is routinely done with classical lamination theory s Classical 
lamination theory uses the four independent material properties of the 
constituent lamina to determine the laminate properties. For a beam, 
the flexural modulus, E, can be cojsputed by E * 12 0^/bh 3 , where h is 
the laminate thickness and b the laminate width. Using the material 
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TAME 6.1 

FINITE ELEMENT INPUT DATA 




Lamination Theory 

Moduli 

Used In 

Dynamic Data For 



Moduli 

Flnite-Clenent 

Finite-Element 





Analysis 

AralysH 



C* 




Lumped 

Initial 

laminate 

Lay-Up 

E c 

£ l 

£ c 

Mass 

velocity 

Number 

(ASl) 


(MSI) 

(Lt») 

(In. /sec.) 

1 

lO!^ 

20.0 

18.5 

18.9 

14.6 

55.6 

207. 

2 

|(15/0-155 5 I s 

19. C$ 

17.22 

18.8 

13.7 



3 

[(30/0/-30) 5 1 s 

15.24 

13.82 

15.9 

8.4 



4 

[(45/0/-45) 5 l, 

11.31 

10.31 

10.9 

5.5 

32.3 

175. 

5 

[(60/0/-60) 5 1 s 

8.88 

8*13 

5.2 

7.1 



6 

({75/C/-7S) 5 ! s 



7.1 

5.9 



7 

[{90/C/-90) 5 1 s 

7.95 

7.27 

7.2 

6.85 

32.3 

175 . 

6 

( (0/90)g| s 

12.06 

10.35 

12.06 

:c.95 



* 

((Cg/Wgij 

18.28 

16.52 

IS. 28 

.6.52 



ID 

UOtyOglj 

£.05 

3.7? 

4,05 

3 . 7 7 



u 

|(4S/-<5/0/90) 4 ) s 

8.91 

8.15 

8.10 

5.56 



12 

l (C/«5/C-4S) 3 /90/O/O l/2 ; s 

14.34 

9.18 

14.34 

9.18 




Material Proaerties Jsed In Uginatg Analysis 


t\ - 2C. 5 MSI ET * 19.5 ms: 

i A 

l\ = l.fc’ MSI EJ * 1.54 MSI 

G « 0.87 MSI 


V 


12 “ 


0.30 
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property data of table 6.1, the oppressive and tensile flexural moduli 
were computed for each laminate type using classical lamination 
theory. These values are presented in table 6.1 and are labeled 
"lamination theory moduli'*. In addition, the flexural moduli of 
laminates 1 through 7 and 11 were empirically determined in ch. 4. The 
values of the moduli at the maxiaiuT curvature levels c? each lar^rate 
type were compared with the lamination theory moduli. The comparisons 
showed poor agreement, particularly the compressive moduli. It was felt 
that the empirically determined moduli were sere accurate. Thus the 


empirically determined moduli were used In the material model for the 
laminates. However, the moduli of laminates 8, 9, 10, and 12 were not 
determined empirically and the lamination theory ssduli were used in the 
material model for those laminates. The flexural moduli used In the 
material model for each laminate type is presented in table 6.1 under 
the column "moduli used in finite elensnt analyses". It is important to 
note that even though the assumptions of classical lamination theory 
were not violated, the lamination theory was pocr at predicting the 


flexural moduli of the laminates. Therefore based on the evidence seen 
here, it appears classical lamination theory cannot be used with 
confidence for the prediction of the large deformation response of the 
laminated beams. 

As stated in ch. 4, the uncertainty in the moduli at the high 
curvature levels was 0.5-1 MSI. The effect of this uncertainty on the 
response is examined in Appendix I. 
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6.4 Finite Element Mesh 

Twenty elements were used to model the 20 in. unsupported portion 
of the specimen between the hinges. The element had a cub^c transverse 
displacement field and a linear axial displacement field. The required 
bean thickness and widths were obtained from table 2. The total load 
eccentricity was equal to the 5/9 in. off-set introduced by the test 
fixture, plus the amount of measured beam camber from table 2. A single 
rigid element was used to model eacn hinge. Figure 6.1 snows the finite 
element mesh used. 


6.5 Static Analysis 

For the static analysis a vertical force vas applied at node 1. 
Ten load steps, from 0 to the maximum load encountered In the 
experiments, were used to compute the static response of the beam. 


6.6 Pynamic Analysis 

To simulate the dynamic load, a lumped mass equal to the mass car 
plus the slider was assigned to node 1. Node 1 n-as giver an initial 
velocity computed from the conservation of momentum of the mass car 
(which had an initial velocity of 235 in. /sec.) Impacting the slider. 
In addition, a constant force In the axial direct -on equal to the weight 
of the lumped mass was applied at node 1. Values of the lumped mass and 
initial velocities used for the three dynamic analyses are given in 
table 6. The starting time step size was 25 u sec. Subsequent tine 
steps were chosen by the program. 
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Figure 6.1 Finite Element Mode! 
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6.7 Static Analysis Results 

The end load -displacement relation from the finite-element analyses 


are displayed with the experimental load-c! ispiacsssent relation '<• fig. 


3.8 for laminate 1, In fig. 3.3 for laminate 2, In fig. 3.11 for 
laminate 3, in fig. 3.14 for laminate 4, In fig. 3.17 for laminate 5, In 
fig. 3.20 fo r laminate 7, In fig. 3.23 for laiJiinate 8, in fly. 3.26 for 
laminate 9. in fig. 3.29 for laminate 10, in fig. 3.32 for laminate 11 
and in fig. 3.35 for laminate 12. For laminates 1 through 8 and 12, the 
finite element analysis agrees well with the experiment up to initial 
failure. However, notice from table 6 the significant difference In the 
flexural moduli comouted from laminate analysis compared with the 
empirical flexural moduli used in the analysis. Fcr 
the !(e/0/-3) 5 l s laminates, if the theoretical flexural moduli had been 
used in the finite element analysis, there would not have been good 
agreement with the experimental load-displacement relation. The finite- 
element program wojld have under predicted the deflections for a given 
load. 


Despite using theoretical moduli, the finite-element comparison for 
laminate type 8 was good. Unlike the previous comparisons, using the 
theoretical moduli in laminate 9 results in an overp-ediction of the 
deflections. For laminate il the finite element analysis underpredicts 
the deflections of the static specimens* ever though the empirically- 
derived flexu**al moduli were used in the analysis. 


6.8 Dynamic Analysis Results 

dynamic analyses was performed for laminates 1, 4, and 7. 
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in the experiments. The initial load spike predicted by an*!**** 


Comparisons of the analyses and the experimental data were made for 
three relations: load-displacement, load-time and displacement-time. 
Analyses were performed up tr the initial failure point in the bear.. 

figure 6.2 shows the comparison of the analytical and experimental 
loac-time relation for laminate 7. The relations are qualitatively 
similar out there are some obvious quantitative differences. In the 
first 0.01 sec. of the load response, much higher frequency components 
of structural vibration are predicted by the analysis than were recorded 

f ho ana 1 % 4 c $ C 

W«'W UIIU ' J J ' J 

had an amplitude of 9845 lbs, more than eight times that which was 
recorded in the experiments. (The load response from the analysis was 
clipped at 4000 lbs. for display in Fig. 6.2), the amplitude of the 
third mode vibratory response was approximately twice that recorded in 
the experiments. However, the frequency of the third mode response was 
predicted quite accurately. Figure 6.3 shews the displacement tise 
relations for laminate 7. Correlation betweer the experiment and 
analysis is good. When the 'oad-time relation and the displacement- time 
relation are cress plotted in figure 6.4 to show the lead-displacesent 
relation for laminate 7, the oscillations in the analysis become 
slightly out of phase with the oscillations in the experiments. 

The same three relations are shown in fig. 6.5 th^j fig. 6.7 f or 


laminate 4 and in fig. 6.8 thru fig. 6. 1C for laminate 1. Comparisons 
similar to those above can be drawn for the correlation of the analysis 
and the experiment for laminate 4 and laminate 1. 

Figure 6.11 shows 5 deformed meshes from the dynamic analysis of 
laminate 4. These deformed meshes were selected for comparison with the 
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Figure 6.2 Dynamic Load vs. T ime: Experimental and Finite 

Element Results for a i <90/0/90) 5 ) s Laitinate 



END DISPLACEMENT, IN. 
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Figure 6.5 Dynamic Load vs. Time: Experimental and Finite 
Element Results for a [ (45/07-45)5)5 Laminate 
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Figure 6.6 Dynamic Displacement vs. Time; Experimental and Finite 
Element Results for a [ (45/0/-45) 5 l $ Laminate 
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Figure 6.7 Dynamic Load vs. Displacement: Experimental and Finite 
Element Results for a ! (45/0/-45) j; 1 * Laminate 
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Figure 6.8 Dynamic Load vs. Time: Experimental and Finite 
Element Results for a 10 1 jq Laminate 
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Figure 6.9 


Dynamic OisDlacement vs. Time: Experimental and Finite 

Element Results for a (0]-> n Laminate 
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Figure 6.10 Dynamic Load vs. Displacement: Experimental and Finite 

Element Results for a (Qijg Laminate 
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5 frames of film shown in fig. 3.4. The time after impact is indicated 
with each frame. The spatial shapes of the deformed meshes in fig. 6.11 
agree quite well with the spatial shapes shown in the film. The times 
after impact are different in figs. 3.4 and 6.11 because there was no 


control in the time increment In the finite element analysis. 

in summary, the analysis qualitatively predicts tne nonlinear large 
deformation resonse of the dynamic tests but quantitatively cverpredicts 


the amplitude of the initial load spike and the subsequent vibration. 
In addition, the analysis predicts mucn higher frequencies in the 
initial 0.01 sec. of the dynamic event. The discrepancies between the 
analyses and the experiments might be due to the lack of structural 
damping in the finite element model. In the experiment there is some 
material damping in the composite specimen, and there is dairp^ng of the 
high frequency components of the load as it is transmitted from the 
slider, to the load transducer, tc the hinge, and finally to the 
composite specimen. The analysis assumes these lead transfers are 
perfectly linear elastic. Also, recall from. ch. ? that the analog 
signal from the force transducer was filtered at 1000 nz before it was 
recorded. To make reasonable comparisons between the experimental data 
and the analysis, the analytical results should have been filtered at 
the same 1000 hz. Filtering analytical output to max® comparisons with 
experimental data is a practice which has been reported oy others 
researching crash behavior (ref. 5). 



Chapter 7 

CONCLUSIONS AND RECOMMENOATICNS 


This study investigated the large deformation bending response of 
simple rectangular cross section composite beams. The study was seen as 
the first step in analyzing the crashworthiness of a composite fuselage 
structure. The overall goals of the study were to: 

1. Design a simple test fixture to Introduce crack-related bending 
loads In beans. 

2. Determine any difference between large deformation static and large 
deformation dynamic response. Specific interest is in the failure 
mode. 

3. Determine the influence of laminate stacking arrangements on the 
dynamic response and failure mode. 

4. Predict the static and dynamic response, using an existing finite 
element program. 


Conclusions 


From this study the following conclusions .may be drawn: 

1. The eccentrically leaded column test fixture proved to be 
successful in introducing large deformation dynamic bending 
loads in to structural elements. The inst-umenUticn recorded 
the deformation and failure response of the beam specimens 
with a minimum of problems. The configuration of the loading 
fixture was such that, under dynamic loading, a high amplitude 
initial load spike and a third mode vibratory response was 
excited in the beam. This more severe dynamic ioac 
environment is desirable for determining the dynamic load- 
deflection behavior of specific structural elements. Such 
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crash conditions will Indeed introduce a severe load 
environment. 

2 . Laminate 9, a {(Qg/90g)!£ laminate and liuiinits 10, a 
[(90 e /0 e )l s laminate exhibited unusual and oistinct failure 
inodes under both static and dynamic loadings. However because 
these laminate had clustered l8Sii"3£, they are not of 
practical interest. 

3. The rest o f the laminate types tested were of some practical 
interest. These laminate types exhibited no differences 
between the failure modes under static and dynamic loading. 

4. The displacement levels at failure and, re’atedly, the strains 
levels at failure were greater for the static tests than the 
dynamic tests. Strain-rate related material properties were 
not thought to be the cause of this result. Instead, the more 
severe load environment experienced by the dynamically tested 


specimens was felt to have beer, the cause. 

5. The progression of failure, l.e., the number cf distinct 

failure events during the 16 in. axial end displacement cf the 
beam, varied with laminate type. Laminate 2, the 
f (15/0/- 15) ^ ] s laminate, had from 3 to 6 failure events with 
associated small drops in the end lead. Laminate 3, the 

( ( 30/0/- 30) $ ! laminate, had from 2 to 3 failure events with 
moderate drops in the load. All other laminates had 1 major 
failure event, producing a drop in the lead 3 f U3 to 305c . 

6. Neither phenomenological (strain tensor polynomial) nor 
mechanistic (maximum fiber strain) failure criteria as 
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successful in predicting the onset of failure. Failure was 
defined as the initiation of a reduction of load carrying 
capability. The strain tensor polynomial is designed to 
predict the onset of any failure such as matrix crack. The 
initiation of a matrix crack did not always coincide with the 
reduction cf load carrying capability. TPs fibers failed ir 
the presence of a matrix crack in an adjacent lamina. The 
matrix crack caused a strain concentration in the fibers and 


8 . 


it is felt this Initiated failure. Since it does not 
accurately measure the fiber strain in the region of failure, 
the strain measured by a strain guage cannot successfully be 
used to predict failure. 

Because composite materials are heterogenous, when a crack 
initiates under bending loads in a lamina, it seldom 
propagates completely through the thickness of the laminate. 
Instead, the crack often turns at a lamina interface, and 
causes del ami nations. This leaves part of the lamina 
undamaged and able to carry a portion of the origins i bending 
load. Failure critera do not address this Issue of the 
prediction of the extent of failure. 

The compressive side surface strain at tbs center of the besn 
was always greater than the tensile side surface strain, even 
tnougn the beam was nearly in a state of pure bending. 
Further investigation revealed that the laminates bad 
blmodular flexural properties. For laminate 4, the 
, (45/0/-45)gJ s laminate, the compressive flexural modulus was 
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measured to be nearly 50% less than the tensile flexural 
modulus. For the unidirectional laminate tne compressive 
flexural modulus was measured to be 23* less than the tensile 
flexural modulus. Because of this, classical lamination 
theory was rot successful in accurately predicting the 
flexural moduli cf the laminates. 

with empirica'ly determined irater'al properties used in the 
material model, the finite element analysis predicted with 
reasonable accuracy the static load-deflection relation, and 
the dynamic load-time, displacement-time and Icad-displicesient 
relations. 


Recommendations 

The state of the art in finite element analysis is such that the 
dynamic structural response can be predicted accurately. However, tne 
program must have available an accurate rated a! aipdel rfhich 
characterizes the flexural rigidity, failure initiation, and failure 
extent. Unfortunately the state of the art in the mechanics of 

- 4 ^ material f rap nn H p*rm< i zip r a omi A 1 HWh/lO 1 Thltt finl x/ 

COfltpOSItc 7law6ria*S Can iUc piw« .gi «-» j **••• j 

empirically determined Infomrafon can be used with ccnf 'cence. 

Therefore further research ir the area of crashworthiness of 
composite structures should be directed J n t.*.c sain areas: 

1. Basic research in the mechanics of composite materials needs 
to be conducted tc being aole to predict the material 
response. The bimodu 1 ar behavior of s lapirate, the 
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initiation of failure in a laminate, and the extent of failure 
once it has initiated need to be studied. 

Further research should be conducted toward empirically 
characterizing the load-deflection behavior of practical 
structural eleirents, such as hat stiffeners or sections of a 
stiffened panel, under both static and dynaiJilc leads. The 
test procedures used in this study would be suitable for such 
a characterization. 
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Appendix A 

EFFECT OF VARIATION IN THICKNESS MEASUREMENT ON THE 
PREDICTED END LOAD-DISPLACEMENT RELATION 

During the cure cycle of the laminated composite plates, excess 
resin from the pre-preg tape flows out of the laminate. The flow of 
resin is not uniform throughout the plate. This results in a variation 
of plate thickness. When specimens are cut from the plate, the 
specimens are of various thickness. The thickness measurements for the 
specimens used In this study were presented in table 2.2. The variation 
in thickness from the end of the specimens to the center was fairly 
consistent and was accounted for in the analytical model. However, 
there was a significant variation in the thickness at the center of each 
oeam, from specimen to specimen, and a random variation in the thickness 
across the width of any given specimen. The uncertainty in the 
thickness measurement was taken to be t 0.002 in. This was one of the 
factors which could account for the deviation ir. the load-displacement 
relation from specimen to specimen. 

To determine the effect of the variation in tnickness on the load- 
displacement relation, the finite element analysis discussed in ch. 6 
was employed. Figure A.l shows the load displacement relation for 3 
unidirectional beams. The center cun/e is from the beam with the 
average center thickness of t = 0.168 in. taken * r om table 2.2. th? ton 
and bottom curve represent the beams with the maximum variation from the 
average center tnickness, t = C.166 in. and t = .170 in. respective’y. 
With the uncertainty of ±0.002 in. in the thickness the^e is an 
uncertainty in the displacement at 3C'G lbs of end load, cf ±0.94 in. So 
with a variation in thickness of less than 1 the variation in 
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Figure A. 1 Effect of Beam Thickness on Stat'c Load-Disc 'acement 
^elation 
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displacement was about 8X. Such a variation in thickness could account 
for scatter of the data from the three replicate specimens. 


Appendix B 

EFFECT OF ECCENTRICITY ON THE LOAD-DISPLACEMENT RELATION 


In the manual construction of a laminated plate, the orientation of 
the lamina can be placed with an accuracy cf, at best, t? . Cess than 
perfect placement of the lamina results in a slightly unsymmetric 
laminate. After the cure cycle the laminate becomes warped by the 
resulting unsymetric residual thermal stresses. The »*rp cf the 
specimens contributes to the eccentricity of the load. The amount of 
warp In the beams of each laminae type was measured as the camber, or 
deviation from perfect straightness, with an uncertainty tC.Ql in. The 
results of these measurements were presented in table 2 . 2 . The amount 
of camber in each laminate type was included in the analytical model. 
However the uncertainty in the camber measurement, taker to be tC.Ol 
in., was one of the factors effecting the deviation in the load- 
displacement relation from specimen to specimen. 

To determine the effect of the minor variations In the eccentricity 
on the load -displacement relation, the finite element analysis discussed 
Ch. 6 was employed. The eccentricity due to the camber of the beams was 
added to the eccentricity from tfis loading configuration sshich -fas 5/8 
in. Figure B.l shows the load-displacement relation for three 
unidirectional beams. The one beam was the average camber, from table 
2.2, of 0.C3 in. Tne other two curves represent the beaas with the 
maximum variation from the average camber, 0.02 and 0.04 in., respec- 
tively. As can oe seen in the figure the three curves are essentially 
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Figure 3.1 Effect o* Beam Eccertr s c i ty or* Static ISfiS^CiSP : 
Relation 
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coincident. Therefore minor variations (*IX) In the eccentricity of the 


load do not significantly effect the lead displaCsnsert rslutio- 


Appendix C 

DIGITAL FILTER ROUTINE 


The random noise in the displacement signal, which was 
subsequently digitized at 4000 samples per second, was smoothed by a 
digital filter routine. A low-pass filter with an ideal cut-off 
frequency of 500 hz. was designed and then programmed in FORTRAN. To 
design the filter, a fast Fourier transform routine transformed a filter 
gain function in the frequency domain to a smoothing function in the 
time domain (ref. 10). Then a window routine computed an optimal 
weighting sequence from the smoothing function. The weighting sequence 
used was 31 clock periods long, i.e., 31 sample points long. 

Table C.l shows the positive portion of the weighting sequence 
used. The table gives the gain characteristics, in decibels, and the 
phase characteristics, in degrees, as a function of frequency. The 

negative portion, I.e., clock period -1 to -15, is the same as the 
positive portion. Thus the filter is synsnetric. 
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TABLE C 

FILTER CHARACTERISTICS 


Clock 

Weight 

Freq 

Sain 

Phase 

Period 

Sequence 

(HZ) 

dbs 

Shift 





(Deg) 


0.12500 

0 

- 0.00000 

0.00000 

0.12142 

250 

-0.00001 

-0.00003 

0.11109 

500 

-6.02061 

-0.00001 

0.09518 

750 

-120.00000 

0.000CG 

0.07544 

1000 

-120.00000 

0.00000 

0.05401 

1250 

-120.00000 

o.ooooc 

0.03307 

1500 

-120.00000 

0.00000 

0.01457 

1750 

-120.00000 

o.ooccc 

0.00000 

2000 

-120.00000 

0.00000 

- 9. 00931 

2250 

-120.00000 

0.00000 

-0.01476 

25CC 

-120.00000 

0.00000 

-0.01543 

2750 

-120.0000C 

o.ooooc 

-0.01294 

3000 

-120.00000 

0.00000 

-0.00876 

3250 

-120.00000 

0.00000 

-0.00440 

3500 

-120.00000 

0.00000 

-0.00118 

3750 

- 120.00000 

o.ooooc 


Appendix D 

STATIC LOAD-DISPLACEMENT DATA 


In this appendix the static load-displacement relation for each 


specimen tested is displayed. 


TL« C i in t MMn rWm- ♦ ^ A 1 n a ^ 
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unloading load-displacement relations. Specimen 10 shows unusual 


behavior because it began to crack immediately upon loading. 
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Figure 0.1 


Static Load-Displacement Relation for Specimens 
1.1 and 1.2 
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Figure 0.2 


Static Load-Displacement Relation for Specimens 
1.3 and 2.1 
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320- « 1S/0/.1S1,], 


240 



ENO OISPLfiCEMENT . IN. 


figure 0.3 Static Load-Displacement Relation for Specimens 
2.3 and 3.1 



END LOAD. LBF END LOAD. LBF 
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Figure 0.4 


Static Load-Displacement Relation for Specimens 
3.2 and 3.3 
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Figure 0.7 


Static load-01 sal acement Relation for Specimens 
5.2 and 5.3 


EM) LOW). LBF EM) LOflO. LBF 
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Figure D.8 


Static Load-Displacement Relation for Specimens 
7.1 and 7.2 


END LORD. LBF ETC LOflO. LBF 
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Figure 0.9 


Static i.oad-0isp1acement delation for Specimens 
7.3 and 8.1 


END LORD. LBF END LORD. LBF 


173 



figure D.1G 


Static Load-DispiacCTent ^elation for Specimens 
8.2 anc 9.3 




END LOflO. LBF END LOflO. LBF 
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Figure D.ll 


Static Load-Displacement Relation for Specimens 
9.1 and 9.2 
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Figure 0.12 


Starfc Load-0 i$p<ac«nent Station for Specimens 
9.3 and 10.1 


END LOAD. LBF END LOAD. LBF 
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O 4 6 12 16 20 

END DISPLACEMENT. IN. 


Figure C.13 


Static wOaa-3isoiace«ent Relation For Specimens 
10.2 and 10.3 


END LOGO. LBF END LOGO. LBF 
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Figure 0.14 


Static load-0 i so lacement Relation for Specimens 
11. 1 and 1.2 







END LOAD. LBF LOAD. LBF 


179 



Figure 0.16 


Static Load-01 sol acement Relation for Specimens 
12.2 and 12.3 



Appendix E 

DYNAMIC LOAD-TIME DATA 


This appendix presents the dynamic load-tine relation for all the 
aminate types tested, with tne excepticr cf laninate type 2 . 
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TIME, SECONDS 

Figure E.i Dynamic Load vs. Time For the 10 ] jq 
l(30/0/-30) 5 j s Laminates 




ENO LORO. LBF 



Figure 1.2 Dynamic Load vs. Time for the ( (45/0/-45)cl 
and ( (60/0/-60)c|_ Laminates 3 
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Figure E.3 Oynawic Load vs. Time for the l(90/0/90)e] s and 
( (0/90) gl s Laminates 







END LOflO. LBF 
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Figure E.5 Dynamic Load vs. Time for the [ (45/-45/0/90)4L and 
[(0/45/0/-45) 3 90/0/0 1/2 ) s Laminates 


Appendix F 

ULTIMATE FAILURE STRAINS 


This appendix presents the ultimate longitudinal strains for al t nf 
the static and dynamic specimens tested. The results are presented In 
Table F.l. The compressive side strain and the tensile side strains are 
recorded in the table. 
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TABLE F.l 
FAILURE STRAINS 


Laminate Specimen c?,,*. t? n . 

Number u,t U1t 



1.1 

.0169 

0.189 

Static 

1.2 

.0153 

.0174 


1.3 

(not strain 

guaged) 


1.4 

.0153 

.0168 

Dynamic 

1.5 

.0146 

.0171 


1.6 

.0152 

.0170 


2.1 

.0086 

.0101 

Static 

2.2 

.0084 

.0096 


2.3 

(not strain 

guaged) 


2.4 

.0075 

.0088 

Oynaailc 

2.5 

.0075 

.0088 


2.6 

(not tested) 


3.1 

.0100 

.0140 

Static 

3.2 

.0104 

.0135 


3.3 

.0104 

.0135 


3.4 

.0091 

.0128 

Dynamic 

3.5 

.0093 

.0127 


3.6 

(not tested) 


4.1 

.0110 

.0150 

Static 

4.2 

.0110 

.0152 


4.3 

.0112 

.0155 


4.4 

.0099 

.0140 

Dynamic 

4.5 

.0100 

.0140 


4.6 

.0100 

.0141 


5.1 

.0122 

.0146 

Static 

5.2 

.0111 

.0134 


5.3 

.0121 

.0148 


5.4 

.0098 

.0118 

Oynamic 

5.5 

.0097 

.0116 


5.6 

.0103 

.0124 
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TABLE F.l (continued) 


Laminate 


Specimen 

e Ult 

c 

e Ult 

Number 




7 


7.1 

(instrumentation problems) 


Static 

7.2 

.0134 

.0140 


7.3 

(instrumentation problems) 



7.4 

.012 

0.130 


Dynamic 

7.5 

.0125 

.0135 


7.6 

.012 

.0135 

8 


8.1 

.0154 

.0170 


Static 

8.2 

.0156 

.0174 



8.3 

.0162 

.0177 



8.4 

.0141 

.0156 


Oynamic 

8.5 

.0141 

.0156 


8.6 

.0150 

.0168 

9 


9.1 

.0137 

.0170 


Static 

9.2 

.0116 

.0132 



9.3 

.0144 

.0176 



9.4 




Dynamic 

9.5 

(no meaningful 

data) 


9.6 



10 


10.1 




Static 

10.2 

(no meaningful data) 



10.3 





10.4 




Dynami c 

10.4 

(no meaningful 

data) 


10.6 



11 


11.1 

(instrumentation problems) 


Static 

11.2 

(instrumentation problems) 



11.3 

.0138 

.0170 



11.4 

.0120 

.0150 


Dynamic 

11.5 

.0123 

.0153 


11.6 

.0120 

.0150 
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TABLE F.i (continued) 


t C 

Laminate Specimen 

Humber 



12.1 

.0127 

.0160 

Static 

12.2 

.0142 

.0176 


12.3 

.0132 

.0161 


12.4 

.0123 

.0160 

Oynamlc 

12.5 

.0135 

.0174 


12.6 

(not tested) 




Appendix 0 

FLEXURAL MODULI DATA 


This appendix presents the data used to calculate the bending 
tensile and compressive moduli of the laminates. Tabulated are the 
axial displacement of the end of the beam, the lateral deflection of the 
center of the beam, the axial load, and the tensile and compressive 
strains. To compute the distance D used in the computation of the 
moduli, 0.725 in. must be subtracted frees the lateral displacement. 
Then the distance from the compressive surface to the neutral surface 
must be added to this difference. (The 0.725 in. represents the offset 
of the vertical reference bar). 
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Table 6.1 

FLEXURAL DATA FOR LAMINATE 1 
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TABLE 6.2 

FLEXURAL OATA FOR LAHINATE 2 


Axial 

Disp. 

(in.) 

Lateral 

Disp. 

(in.) 

i o ■ mIIIII 

Tensile 
Strain (e ) 
M 

Compessive 
Strain (e 
(we) 

0.50 

2.73 

129 

1820. 

2020. 

1.00 

3.60 

153 

2940. 

3350. 

1.50 

4.25 

166 

3830. 

4430. 

2.00 

4.76 

175 

4580. 

5350. 

2.50 

5.25 

181 

5270. 

6170. 

3.00 

5.65 

187 

5900. 

6920. 

3.50 

6.02 

192 

6500. 

7630. 
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TABLE 6*3 

FIEX’JRAL DATA FOR LAMINATE 3 
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Table G.4 

FLEXURAL DATA FOR LAMINATE 4 


Axial 

Disp. 

(in.) 

Lateral 

Disp. 

(In.) 

Axial 

Lead 

(lbF) 

Tensile t 
Strain (e l ) 
(ve) 

Compressive 
Strain (e 
(we) 

0.50 

2.75 

57 

1850. 

1950. 

1.00 

3.60 

67 

2940. 

3270. 

1.50 

4.25 

73 

3790. 

4360. 

2.00 

4.80 

76 

4490. 

5310. 

2.50 

5.25 

79 

5100. 

6200. 

3.00 

5.69 

82 

5650. 

7020. 

3.50 

6.04 

84 

6160. 

7780. 

4.00 

6.40 

36 

6640. 

8520. 

4.50 

6.70 

88 

7090. 

9230. 

5.00 

6.99 

90 

7530. 

9900. 

5.50 

7.25 

92 

7940. 

10560. 

6. CO 

7.50 

94 

8340. 

11200. 

5.50 

7.71 

96 

8730. 

11830. 

7.00 

7.95 

98 

9110. 

12440. 

7.50 

8.15 

ICO 

9480. 

13050. 

8.00 

3.30 

1C2 

9850. 

13650. 

8.50 

8.50 

103 

10200. 

14250. 

9.00 

8.67 

105 

10500. 

14830. 



TABLE $.5 

FLEXURAL OATA FOR LAMINATE 5 
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TABLE G.6 

FLEXURAL DATA FOR LAMINATE 6 


Axial 

Disp. 

(In.) 

Lateral 

Disp. 

(in.) 

Axial 

Load 

(lbF) 

Tensile* 
Strain (e ) 

(ue) 

Compressive 
Strain (e c 
(-0 

0.50 

2.75 

53 

1770. 

1930. 

1.00 

3.60 

62 

2920. 

3200. 

1.50 

4.25 

66 

3830. 

4220. 

2.00 

4.80 

69 

4610. 

5110. 

2.50 

5.25 

70 

5310. 

5900. 

3.00 

5.68 

72 

5960. 

6660. 

3.50 

6.03 

74 

6540. 

7340. 

4.00 

6.38 

78 

O 

o 

8000. 

4.50 

6.70 

77 

7640. 

8640. 

5.00 

6.98 

78 

8140. 

9260. 

5.50 

7.25 

80 

3710. 

9880. 

6.00 

7.50 

ac 

9340. 

10490. 

6.50 

- • V 

31 

9820. 

11060. 

7.00 

7.95 

32 

10250. 

1159C. 

7.50 

8.15 

83 

10670. 

12140. 

8.00 

8.35 

84 

11100. 

12360. 
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TABLE G.7 

FLEXURAL DATA FOR LAMINATE 7 


Axial 

Oisp. 

(in.) 

Lateral 

Disp. 

(in.) 

Axial 

Load 

(lbF) 

Tensile ¥ 
Strain (t z ) 
M 

Compressive- 
Strain (c C 
(ue) 

0.50 

2.75 

53 

1800. 

1920. 

1.00 

3.60 

62 

2930. 

3110. 

1.50 

4.25 

67 

3830. 

4040. 

2.00 

4.80 

71 

4610. 

4840. 

2.50 

5.25 

73 

5310. 

5560. 

3.00 

5.68 

76 

5960. 

6220. 

3.50 

6.03 

78 

6560. 

6840. 

4.00 

6.38 

79 

7130. 

7440. 

4.50 

6.70 

81 

7680. 

8010. 

5.00 

6.98 

83 

8190. 

6560. 

5.50 

7.25 

85 

8690. 

9090. 

6.00 

7.50 

86 

9180. 

9610. 

6.50 

7.75 

88 

9680. 

10110. 

7.00 

7.55 

89 

10130. 

1C610. 

7.50 

8.15 

91 

10530. 

11110. 

8.00 

8.35 

92 

10970. 

11580. 

8.50 

8.50 

94 

11370. 

12060. 

9.00 

8.68 

95 

11830. 

12550. 
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TABLE 6.8 

FLFXURAL DATA FOR* LAMINATE II 


Axial 

Oisp. 

(In.) 

Lateral 

Dlsp. 

(in.) 

Axial 

Load 

(lbF) 

Tensile f 
Strain (e l ) 
(-«> 

Compressive 
Strain (e 1 
(ue) 

0.50 

2.72 

72 

2080. 

2140. 

1.00 

3.55 

83 

3320. 

3480. 

1.50 

4.22 

88 

4310. 

4580. 

2.00 

4.78 

92 

5150. 

5550. 

2.50 

5.20 

95 

5890. 

6430. 

3.00 

5.65 

98 

6560. 

7250. 

3.50 

6.01 

100 

7180. 

8020. 

4.00 

6.35 

102 

7760. 

8750. 

4.50 

6.70 

104 

8310. 

9460. 

5.00 

6.98 

106 

8830. 

101 4C. 

5.50 

7.25 

107 

9330. 

10810. 

6.00 

7.50 

109 

9810. 

11450. 

6.50 

7.7? 

111 

1029G. 

12 ICO. 

7.00 

7.95 

112 

10730. 

12700. 

7.50 

8.14 

114 

11170. 

13310. 

8.00 

8.31 

116 

11590. 

13920. 

8.50 

8.50 

118 

12020. 

14520. 

9.00 

8.66 

119 

12450. 

15110. 





Appendix H 


EFFECT OF TRACKING THE LOCATION OF THE NEUTRAL SURFACE 
ON THE LOAO-O ISPLACEMENT RESPONSE 


As described in ch. 4, the laminates exhibited bimodular 
properties. With diffe-ent tensile and compressive flexural moduli, the 
neutral surface of the beam was no longer coincident with the 
midplane. For the analysis performed in ch. 6, the neutral surface and 
the midplane were assured to be coincident. "he effect of that 
assumption on the static load-displacement response is examined here. 

To examine this effect, laminates I and 3 were studied. Laminate 3 


showed the greatest bimodu larHy. Tne compressive flexural modulus was 
50% less than the tensile flexural modulus, and hence, this laminate 
should shew the largest effect of the shifting of the neutral surface. 
Laminate I showed a more moderate bi modularity with the compressive 
moduli being 23 % less than the tensile moduli. Figures H.l and H.2 show 
the load-displacement relation for laminates 3 and i, respectively. In 
the upper curve the analysis assumes the neut^a 1 surface to be at the 
midplane. In the lower curve the analysis tracks the location of the 
neutral surface. For laminate 3 the error is significant. At 150 lb of 
force, there is a difference of 2 in. in the predicted axial 
deflection. This is a 25% error. For laminate L the error is much 
less. At 300 lb of force, there is a difference of 0.25 in. in the 
predicted axial deflection. This is an e’ror of only 2 %. 
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Appendix I 


EFFECT OF UNCERTAINTY IN THE MODULI ON THE 
LOAD-DISPLACEMENT RELATION 


As mentioned in ch. 5 tne uncertainty in the empirical moduli 
calculations was about ± 1 MSI. This is one factor which could account 
for the deviation between the analytical results and the experimental 
data. To examine the effect of the uncertainty in the moduli or the 
load-displacement relation, the finite slsssnt analysis wss ejnployeu. 

Figure 1.1 shows the load-displacement relation for 3 
unidirectional beams. The center curve has the flextral properties 
taken from Table 6.1. The other two curves represent the range of the 
uncertainty in the flexural properties. Uncertainty in the modulus of 
±1 MSI results in a uncertainty in the displacement of ±1.65 in., at the 
12 in. displacement level. """ms represents a 13% dsvisticr »n the 
predicted response. 
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